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ABSTRACT 

We present new proper motions from the 10 m Keck telescopes for a puzzling 
population of massive, young stars located within 3'.'5 (0.14 pc) of the supermas- 
sive black hole at the Galactic Center. Our proper motion measurements have 
uncertainties of only 0.07 mas yr~^ (3 km s~^), which is > 7 times better than 
previous proper motion measurements for these stars, and enables us to measure 
accelerations as low as 0.2 mas yr~^ (7 km s~^ yr"^)- Using these measurements, 
line-of-sight velocities from the literature, and 3D velocities for additional young 
stars in the central parsec, we constrain the true orbit of each individual star and 
directly test the hypothesis that the massive stars reside in two stellar disks as 
has been previously proposed. Analysis of the stellar orbits reveals only one of 
the previously proposed disks of young stars using a method that is capable of 
detecting disks containing at least 7 stars. The detected disk contains 50% of the 
young stars, is inclined by ~ 115° from the plane of the sky, and is oriented at 
a position angle of ~ 100° East of North. Additionally, the on-disk and off-disk 
populations have similar K-band luminosity functions and radial distributions 
that decrease at larger projected radii as oc r~^. The disk has an out-of-the-disk 
velocity dispersion of 28 ± 6 km s~^, which corresponds to a half-opening angle 
of 7° ± 2°, and several candidate disk members have eccentricities greater than 
0.2. Our findings suggest that the young stars may have formed in situ but in a 
more complex geometry than a simple, thin circular disk. 
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Introduction 



The cen ter of our Galaxy har 
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3 ors not only a supermassive black hole (Sgr A* , M, 
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Ghez et al. 1998, 200G, 2003, 2005b 



Eisenhauer et al.l l2005l ) , but also a population of 



massive (10-120 Mq ), young (<10-100 Myr) stars whose existence is a puzzle. The origin of 
such young stars has been difficult to explain since the gas densities observed today are orders 
of magnitude t oo low for a gas clump to overcome the extreme tidal forces an d collapse to 
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form stars (e.g. ISanderslll992l : lMorrislll993l : iGhez et al.ll2005bl : lAlexanderl l2005l . for reviews 



And yet, within the central parsec of our Galaxy, nearly 100 stars have been classified as 
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more evolved massive stars having ages as young as 6±2 Myr ( 


Paumard et al. 


2006). 



ulations of young st ars have also been observed in the nuclei of other galaxies, such as M31 



(iBender et al.ll2005l ). suggesting that star formation near a supermassive black hole may be 
a common, but not understood, phenomenon in galaxy evolution. The close proximity of 
the black hole at the center of the Milky Way provides a unique laboratory for studying 



this " paradox of youth" (e.g. iGhez et al.ll2003l . l2005bl : ISchodel et al.ll2003l : Eisenhauer et al. 



20051). 



Proposed resolutions to the paradox of youth can be grouped into several broad cate- 
gories, including (1) rejuvenation of an older population such that older stars appear young, 
(2) dynamical migration from larger radii, and (3) in situ formation. Rejuvenation sce- 
narios include stripping ( Davies et al. 19981: Davies fc King 2005) or tidal heating of the 
atmospheres of old stars ( Alexander fc MorrisI 20031 ). or combining multiple low mass stars 



via colhsional mergers to form a higher- mass hot star akin to a "blue straggler" (jLedll996 



Morrislll993l : iGenzel et al.ll2003l ). Although these processes may be candidates for explaining 
the closest young stars within the central arcsecond, they cannot account for the OB giants, 
OB supergiants, and Wolf-Rayet stars that are located at larger radii (1"-14"), since the 
rate of collisions is too low to produce the observed total numbers. Thus, it appears that 
these massive young stars must have formed, or were deposited, in the central region within 
the last 4-8 Myr. Dynamical migration scenarios attempt to resolve the paradox of youth 
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with the formation of a massive star cluster at larger distances from the black hole (3-30 
pc). Such a cluster would spiral in d ue to dynamica l friction and deposit stars at smaller 
radii where they are observed today (iGerhardl l200ll ). However, for a cluster to reach the 
central parsec in only a few million years, it must be very massive and centrally concen- 



trated (iKim fc MorrisI 120031 : iPorteries Zwart et al.ll2003l : McMillan fc Porteries ZwartI 12003 
Giirkan fc Rasiolboosh . and it may even req uire the existence of an intermediate-mass blac k 



hole (IMBH) as an anchor in the cluster core (IHansen fc Milosavljevidl2003l : iKim et al.ll2004j ). 
In situ star formation scenarios can resolve the par adox of youth if a massive, s elf-gravitating 
gas disk was once present around the black hole (ILevin fc Beloborodovll2003l ). Such a disk 
would be sufficiently dense to overcome the strong tidal forces, and gravitational instabilities 
would then lead to fragmentation and the formation of stars, as has been suggested in the 
context of b oth the Galactic Center circumnuclear disk and AGN accretion disks in other 
galaxies (e.g. Kolvkhalov fc Svunvaev 198ol: Shlosman fc Begelrnan 1989 : Morris fc Serabyn 
19961 : ISanderj|l998l : lGoodmanll2003l : iNavakshin fc Cuadralboosl l 



Insight into the origins of the massive, young stars may be obtained through obser- 
vations of the spatial distribution and stellar dynamics of this population. Already, high- 
resolution infrared imaging and spectroscopy h ave shown that the young stars between O'.'S 
and 14"(0.02-0.6 pc) exhibit coherent rotation (iGenzel et al.ll2000l ). Analyses of the statis- 
tical properties of the three-dimensional velocity vectors for these stars suggest that they 
may reside in two disks. The fi rst proposed disk has a cloc kwise sense of rotation, as pro- 
jected onto the plane of the sky (ILevin fc Beloborodovll2003l . hereafter: clockwi se-rotating or 



CW disk) , while the second proposed disk is counter-clockwise-rotating (CCW iGenzel et al. 



20031 ) a nd is nearly perpendi cular to the first. The proposed disks extend from ~0'.'8 to at 
least 7^^ (jPaumard et al.ll2006l ). Other velocity vector analyses show that there are possible 
co-moving groups or clusters of st ars, including the IRS 13 cluster, which is proposed to 
lie within the putative CCW disk jMaillard et allbooi ISchodel et al.lboosh . and the IRS 



16SW co-moving group, which are also consistent with the proposed CW disk ( ILu et al. 



20051 ). The two proposed disks are inferred to be oriented with an inclination and angle to 
the ascending node of [icw=127°± 2°, Qcw=99°± 2°] and [iccw=2A°± 4°, Qccw=lQ7°± 
7°] and to have a finite angular thickness of A6cw ~ 14° and A6cw ~ 19° where A6 is 
the standard deviatio n of the orbital inclinations distributed normally about the disk plane 
(jPaumard et al.ll2006l ). The thickness of the stellar disks has been attributed to thickening 
as a result of gra vitational interactions between the two disks, which provides an estimate of 
the disk masses ( iNayakshin et al.ll2006l ). The derived mass is smaller than the mass inferred 
from the nu mber of observed young stars, assuming a Salpeter initial mass function (IMF); 
accordingly, INayakshin et al.l (120061 ) suggest that the disks have a top-heavy mass function. 
Both in situ gas disk and in-spiraling star cluster formation scenarios have been used to 
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explain the kinematics of this young star population and to predict that the stars should lie 
in a common orbital plane. However, the presence of two stellar disks with similarly aged 
populations requires either two nearly concurrent gas disks or two infalling star clusters; and 
both of these scenarios are difficult to produce. Therefore, to understand the recent star 
formation history, it is critical to measure the orbital planes of individual stars in order to 
confirm the existence of the two stellar disks previously derived from a statistical analysis of 
velocity vectors alone. 

The in situ gas disk and inspiraling star cluster formation scenarios predict different 
structures and evolutions for the resulting stellar disk, particularly with respect to the ec- 
centricities and radial distribution of stars within the disk. Early models of a self-gravitating 
gas disk around the supermassive black hole at the center of the Milky Way produce stars 
with a steep radi al profile in the disk surface density, S oc r°, with a ~ —2 (ILin fc Pringle 
19871 : iLevinI 120071 ). These models typically result in stars on circular orbits as would be the 
case for the slow build up of a gas d isk that is circularized before there is sufficient mass for 



gravitationa l instabilities to set in (IMilosavljevic fc Loebl l2004l : iNayakshin fc Cuadral 12005 
Levin! 120071 ) . The stellar eccentricities of an initially circular disk can relax to higher ec- 

0.15 for a normal IMF or e 



centri cities up to a 



rm-t 



rms ~0.3 for a top-heavy 
More recent models have also shown 



IMF flAlexander et al.l 120071 : ICuadra et all 120081 ) 
that star formation can occur rapidly before circularization in an initially eccentric disk 
as mig ht result frorn the infall of a s i ngle r aassive molecular clou d or a cloud-cloud col- 
lision jSandersI 1 19981 : iNavakshin et all I2OO7I : [Alexander et al.l I2OO8I ). These eccentric self- 
gravitating accretion disk models typically produce a more top-heavy IMF than initially 
circular disks. On the other hand, an inspir aling star cluster would d issolve into a disk of 



stars with a fiatter radial profile (S oc r iBerukoff fc Hansen! |2006| ) whose orbital eccen 



or eccentric ( 


Porteeies Zwart et al. 


2OO3I: 


McMillan & Porteeies Zwart 


2003; 


Kim & Morris 


2003; 


Kim et al. 


2004; 


Giirkan & Rasio 


2005; 


Berukoff & Hansen 


2006 


). Previous measure- 



ments of the radi al distribution of you ng stars yields a steep radial profile consistent with 
in situ formation (jPaumard et al.ll2006l ). Also, the eccentricities of the stars have previously 
been estimated from observations by assuming that the stars orbit in a disk; however, there 
are confiicting results claiming that the stars in the cloc kwise-rotating disk are o n nearly 
circular orbits (jPaumard et al.ll2006l ) or on eccentric orbits (IBeloborodov et al.ll2006l ) . Deter- 
mining the radial profile and stellar eccentricities of stars in a disk may provide observational 
constraints on the origin of the young stars. 

We present an improved proper motion study that yields an order of magnitude more 
precise proper motions and the first measurement of accelerations in the plane of the sky for 
stars outside the central arcsecond. By combining the stellar positions, proper motions, radial 
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velocities, and accelerations, we estimate stellar orbital parameters and test whether the 
young stars reside on one or two stellar disks in a more direct manner than previous methods 
using only velocity information. This provides a direct test of the existence, membership, 
and properties of these disks. The observations are described in ^ and the astrometric 
analysis procedure and results are detailed in ^ Orbit analysis and results are presented in 
^ and ^ and a discussion of the implications for the origin of the massive, young stars at 
the Galactic Center is presented in ^ 



2. Observations 



This study utilizes 29 epochs of high-resolution, infrared images of the Galaxy's central 
stellar cluster, which were taken from 1995 to 2005 using both speckle and laser guide 
star adaptive optics (LGS AO) observing techniques on the W. M. Keck 10 m telescopes. 
These data sets are listed in Table [T] and all but the additional LGS AO observation from 



2005 are described i n detail in earlier papers (iGhez et al.lll998l . l2000l . l2005bl : ILu et al.l 12005 
Rafelski et al.l 120071 ). Columns 3 and 4 list the individual exposure times and the total 



number of frames for each epoch of data. All 2 7 speckle imaging observations were taken 



using the facility near-infrared camera, NIRC (IMatthews fc Soiferl Il994l : iMatthews et al. 



19961 ). which has a plate scale of ~20 mas per pixel, and a 5'.'22 x 5'.'22 field of view. The 
two adaptive optics imaging observations used the facility LGS AO system (jwizinowich et al" 



20061 : Ivan Dam et al.l 120061 ) and the near-infr ared camera, NIR C2 (PI: K. Matthews) with 
a plate scale of 9.963 ± 0.006 mas per pixel Jchez et allboosh and a 10'.'2 x 10'.'2 field of 
view. While the laser guide star is used to correct most of the atmospheric aberrations, the 
low-order, tip-tilt terms were corrected using visible observations of USNO 0600-28577051 
(R = 13.7 mag and ArsgrA* = 19"). 

In addition to the 27 speckle observation and the 2004 LGS AO observations described 
in previous works, a new LGS AO data set was obtained in 2005 June. This data set was 
taken using two different narrow-band filters, Kqo (Ao=2.289 /im, AA=0.027 /im) and Kcont 
(Ao=2.270 ^m, AA=0.030 /im), rather than the K' broadband filter used for the 2004 LGS 
AO observations. For each filter, images were taken in a 5 position pattern around a 4'.'0 
box with exposure times of 36 s (texp = 7.2 s, 5 coadds) and 59.5 s (texp = 11-9 s, 5 coadds) 
for the Kco and Kcont filters, respectively. The choice of narrow-band filters was driven by 
a different project and the data sets from the two filters were combined together for the 
present study (see §3.ip . Resulting Strehl ratios were ~0. 25-0. 35 in the individual frames. 
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3. Astrometric Data Analysis and Results 



The goal of this analysis is to obtain high precision astrometry for a sample of young 
stars that are candidate disk members and have existing radial velocity measurements. Based 
on spectroscopic identif ication, there are curr ently 90 known young stars with radial velocity 
measurements listed in iPaumard et al.l (120061 ) based on high quality ( "quality 1 or 2" ) spec- 
tral classifications. We define a primary sample that includes those known young stars found 
in our astrometric data sets that have projected radii between O'.'S and 3'.'5. The inner radius 
is set by the proposed inner edge of the clockwise disk of young stars and young stars inte 



riorto this radius are on more randomly oriented orbits (IGhez et al.ll2005bl : lEisenhauer et al. 



2005h. The out e r radi us is set by the field of view of the speckle data sets. Over this region, 
Paumard et al.l (120061 ) note that all young stars brighter than K=13.5 should be identified, 
which includes OB giants and supergiants. A total of 32 such young stars are in our 11 
year astrometric data set and comprise the sample for this study. Of the 32 stars in our 
sample, 23 are among the 36 stars thought to be part of the clockwise disk, 2 are among the 
12 candidate members of the count er-clockwise disk, and the remaining 7 are among the 42 
stars not assigned to either disk by lPaumard et al.l (120061 ). 



We also define an extended sample that includes both the primary sample of 32 stars and 
an additional 41 young stars found by iPaumard et al.l (120061 ) at larger radii that are outside 
the field of vi ew of our astrorn e tric m easurements. The astrometry for the additional 41 stars 
is taken from Paumard et al. ( 2006 1^. which has an order of magnitude lower precision and 
lacks any constraints on the accelerations. However, we use the extended sample to explore 
the kinematics of the young stars at larger radii with the same analysis techniques used on 
the primary sample. We also note that the spectroscopic observations used to identify the 
young stars at larger radii were taken in a different setup than in the central regions, with 
lower spectral resolution and lower Strehl; thus the completeness limit may be somewhat 
brighter in this region. However, any difference is statistically insignificant given that a two- 
sample KS-test yields a 50% probability that the primary sample and those additional stars 
added to the extended sample have the same K-band luminosity function. The extended 
sample is used only to supplement our analysis; therefore, to avoid confusion, all analysis 
and results are reported for the primary sample, which has more precise proper motions and 
accelerations, unless specifically noted otherwise. 

Astrometric positions for the young stars in the primary sample are extract ed from the 



imaging data sets listed in Table [T] using similar techniques to those described in lGhez et al. 



^We note that there are 4 additional young stars a t larger radii that are not included in our extended 
sample since they do not have proper motions listed in IPaumard et aL ( 20061 ). 
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(119981 . mm, ILu et alJ (120051 ). and IChez et all (l2005bf ). with the following key changes: 
(1) geometric distortion is corrected in the speckle images using an improved distortion 
solution (see §3.11 Appendix Rl), (2) speckle imag es are combined with an improved algorithm 
developed and implemented by iHornsteinI (120071 ). and (3) image coordinates are transformed 
between data sets with more degrees of freedom (see §3.2p . Sections |3. 11 and [3^21 describe the 
analysis in detail and Section 13.31 presents the astrometric results. 



3.1. Image Processing 



To achieve precise astrometry, the basic image reduction steps, particularly geometric 
distortion correction, must be carefully implemented. First, both speckle and LGS AO 
individual exposures are processed using standard techniques of sky subtraction, fiat-fielding, 
and bad pixel correction. Next, the images are transformed to correct for optical distortion. 
For the LGS A0/NIRC2 images, optical d istortions are well characterized at the ~2 milli- 
arcsecond level over 2" (IGhez et al.l |2008| . Appendix by the pre-ship review distortion 
coefficient and the dist ortions are removed from the images using the IRAF routine. Drizzle 
( iFruchter &: Hookll2002l ). The speckle images, obtained wi t h NIR C, have a known off-axis 
distortion that can be corrected as described in lGhez et al.l (Il998l ). However, this distortion 
solution does not account for any distortion introduced by the additional optics in the NIRC 
reimager, which magnifies the image scale by a factor of ~7 from seeing limited sampling to 
diffraction limited sampling. Speckle data sets were acquired in such a way as to minimize 
the effects of this residual distortion in the center of the field of view and have resulting 
residual distortion errors that are smaller than the typical centroiding error, which is ~2 
mas, for stars at radii < O'.'S. However, astrometric uncertainties for stars outside this region 
are dominated by the uncorrected distortion, which grows to ~6 mas near the field edge at 
a radius of 2'.'5 (IGhez et al.ll2005bl ). In order to characterize the residual distortion in NIRC, 
simultaneous images of the Galactic Center were obtained with both NIRC and NIRC2 with 
the NIRC2 images serving as a reference coordinate system (see Appendix . The speckle 
image distortion is mapped by comparing stars' positions in both NIRC and NIRC2 images. 
As shown in Appendix |Al the resulting NIRC to NIRC2 transformation is characterized at 
the ~2 mas level over the entire field of view. 



After distortion correction, individual exposures are combined into a final diffraction- 
limited image using different methods for speckle and LGS AO data sets. Speckle images 
are produced by first rejecting the low Strehl ratio frames (typically 75% of frames are 



^http: / /www2. keck. hawaii.edu/inst/nirc2/ 
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rejected) and then stack ing the remaining frames using a weighted shift-and-add (SAA) 
routine (IHornsteinl 120071 ). The resulting combined images have a point-spread function 
(PSF) composed of a diffraction-hmited core (FWHM~0'.'055) on top of a broad seeing 
halo (FWHM~0'.'4). The improved image combination algorithm attempts to maximize the 
signal-to-noise ratio (SNR) of the final image while preserving the highest spatial resolution. 
Quantitatively, the weighted SAA method doubles the fraction of light contained in the 
diffraction-limited core (fror n 3.5% to 7.0%) over the standard SAA scheme with no weight- 
ing and no frame rejection (IHornsteinl 120071 ) . The LGS AO individual exposures are all of 
similar quality and are thus all averaged together, without weighting, in order to produce 
the final high-resolution image for each data set. Although the 2005 June data were taken 
in two different filters [Kco and Kcont), all the images were combined together to increase 
the final SNR. While photometry from this epoch is marginally impacted, the astrometry is 
comparable to other epochs. Each data set was also sub-divided to produce three equivalent 
quality (randomized in time) subsets to make three images used for determining photometric 
and astrometric uncertainties. The resulting images are summarized in Table (H including 
the achieved spatial resolution (FWHM) and the Strehl ratio. 



3.2. Stellar Positions and Coordinate Transformations 

In order to extract astrometric information for the sample of young stars, the coordinate 
system from each data set is transformed into a common reference frame using the stars in 
each image to determine the transformation parameters. Since the accuracy of this trans- 
formation relies on the assumption that there is no net rotation of the sample, we use all 
stars detected in each data set, not just the young stars, in this analysis. The steps for (1) 
measuring stars' positions in each epoch, (2) transforming to a common (relative) reference 
frame, and (3) determining the absolute coordinate system are described below and utilize 
all stars detected in the data sets; then as a final step, the young star sample is extracted. 

In each data set, stars are identified and their positions meas ured using the IDL point- 



spread function fitting routine "StarFinder" (jPiolaiti et al.ll2000l ). StarFinder generates a 



PSF from several bright stars in the field and cross-correlates the resulting PSF with the 
image. The PSF was iteratively constructed using IRS 16C and IRS 16NW for the speckle 
maps and IRS 16C, 16NW, 16NE, 16SW, 33E, 33W, 7, 29N, and GEN-f2.33-F4.60 for the 
LGS AO images. Candidate stars are those for which StarFinder correlation peaks have a 
correlation value higher than 0.8 and positions and fluxes are extracted by fitting the PSF to 
each correlation peak. From the candidate star list, spurious detections are then eliminated 
by requiring that each star be detected in all three of the subset-images with a correlation 
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of higher than 0.6. The positional centroiding uncertainties for each candidate star are 
estimated from the rms of their locations in the three subset-images, and an additional 
systematic error term of 0.88 mas is added in quadrature to all stars in LG S AO epochs to 



account for residual distortion in the central 5" of NIRC2 (iGhez et al.ll2008l ). The candidate 
stars are flux calibrated using the apparent magnitudes of the non-variable stars, IRS 16C, 
IR S 16SW-E, S2-17, Sl -23, Sl-3, Sl-4, S2-22, S2-5, Sl-68, SO-13, and Sl-25, as measured 



by lRafelski et al.l (120071 ). The star detections from each epoch are cross-identified with stars 
from all other epochs and those stars that are detected in at least 16 out of 29 epochs are used 
to create a master star list. The threshold of 16 or more epochs is used in order to insure high 
astrometric precision; for a threshold of less than 16 epochs, the number of detected stars 
rises dramatically as does the number of sources showing significant (^Sc) accelerations in 
non-physical directions, indicating a high frequency of false detections (see §3.31 for further 
discussion). Stars in the master list are also examined for source confusion, which may occur 
when two stars pass close enough to each other such that StarFinder only detects a single 
source with biased astrometry rather than detecting both stars. Source measurements from 
individual epochs are rejected if two stars pass within 55 mas (~1 spatial resolution element) 
of each other and only one source is detected by StarFinder. The results of this stage of the 
analysis are summarized in Table [1], which provides for each data set the average centroiding 
error for the brightest stars (K<13; also see Figured]) and the sensitivity as estimated by the 
peak in a histogram of the K-band magnitudes (bins = 0.1 mag) of all the stars in the data 
set. Averaged over all stars in all maps, the centroiding uncertainties have a mean value of 
1.6 mas for the brightest stars (K < 13 mag) and 3.4 mas for the fainter stars (13 < K < 16 
mag). 

The coordinate system for each image is transformed to a common local reference frame 
defined by the 2004 July LGS A0/NIRC2 image's coordinates and pixel scale. This particular 
LGS AO epoch was chosen as the reference because the NIRC speckle distortion solution is 
tied to this epoch, thus providing a smooth transition between speckle and LGS AO data 
sets. The procedure for deriving the coordinate transformation for all of the data sets is 
non-trivial, since the stars in the images have detectable motions. Optimal alignment is 
ac hieved by m i nimizi ng the error-weighted, net displacement for all the stars as described 



by iGhez et al.l (119981 ) while allowing for translation, rotation, and two magnifications in 
arbitrary, but perpendicular, directions. This is a higher order transformation than was 
used in our earlier astrometric works, which only allowed for translation and rotation. The 
new transformation equations have the form 

Xpix = ^0 + CLix'pi^ + a2ypix (1) 
Vpix = bo + hy'pi^ + foa^ix (2) 
where Xp^^ and j/^j^. are the input detector coordinates in pixels and Xpi^ and i/pix are the 
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output coordinates for each star, and all o ther variables are fr ee parameters that are common 
across all stars in the alignment fit. As in lGhez et al.l (l2005bl ). stars within O'.'S of Sgr A* are 
excluded from the transformation as they exhibit large non-linear motions. Additionally, all 
spectroscopically ident ified young stars ar e excluded from the transformation as they have 
a known net rotation (IGenzel et al.ll2000l ). Initially, each image is aligned to the reference 
image by assuming the stars have no proper motions and finding the best-fit values for the 
free parameters of the transformation, ao, Oi, 02, 6o; ^i; ^2, for that image. However, after a 
first pass at the alignment of all the images, proper motions are estimated and used to refine 
the alignment solutions in a second pass. Sources with estimated proper motions higher than 
1.5 mas yr~^ (600 km s~^) are excluded from the transformation resulting in the elimination 
of 2 sources that are near the edge of the speckle field-of-view and suffer from edge effects. 
Alignment uncertainties are estimated by a half-sample bootstrap method (iBabu fc Feigelson 



19961 : iGhez et al.l l2005bl ) and are small (~0.2 mas for stars at r <2") compared to the 



centroiding uncertainties (see Figured]). Alignment and centroiding uncertainties are added 
in quadrature to produce a final relative positional uncertainty for each star at each epoch. 
The resulting astrometric data set contains stellar positions and uncertainties for all epochs, 
transformed into the 2004 July NIRC2 pixel coordinate system {xpix,ypix)- 

The relative positions and uncertainties are transformed into J2000 absolute astrometric 
coordinates defined by radio observations of SiO masers and Sgr A*. Using observations of 
the SiO masers in the infrared, a set of infrared ab solute astrometric standards are defined 
in a process described in detail in an appendix of iGhez et al.l (120081 ). These astrometric 
standards are used to derive the transformation from 2004 July NIRC2 pixel coordinates into 
absolute coordinates. A statistically insignificant adjustment is made to place the origin at 
the dynamical center of SO-2's orbit, which is known to high precision, by offseting from the 
radio position of Sgr A* by 1 mas to the East and 5 mas to the South. This offset is well within 
the absolute astrometric uncertainty of ~6 mas for Sgr A* (IGhez et al.ll2008l ). The stellar 
positions in all epochs are thus expressed in arcseconds offset from the dynamical center with 
+x increasing East and +y increasing North and can be converted into celestial coordinates 
using (x, y) = (cos 5 Aa, A5) §. Positional uncertainties are taken as the quadratic sum of 
the relative errors, which dominate, and the absolute error from uncertainties in the plate 
scale and position angle. Errors in the relative position of Sgr A* (~2 mas) are incorporated 
later during the orbit analysis stage as a parameter of the potential of the supermassive 
black hole (see From the resulting absolute astrometric data set, the sample of young 
stars is extracted. 



■^When converting from (x, y) to (Aa, A6), higher order terms are neghgible (0.06 mas over 5") because 
the celestial sphere is sufficiently flat over our field of view. 
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3.3. Proper Motions and Acceleration Results 

For each of the young stars in the sample, positions, velocities, and accelerations in 
the plane of the sky are derived by fitting second-order polynomials to the star's position 
as a function of time, weighted by the positional uncertainties. The polynomials are fit 
independently in x and y coordinates and have the form 

= Xref + Vx^refii ^re/) "I" '^(^x,ref{t tref) (3) 
y{t) = yref + Vy^refit " Uef) + ^^y^refit " Ueff (4) 

where t is the time in years, tref is a reference time taken to be the mean of the time of all 
epochs weighted by the positional uncertainties for each star, Xref and yref are the positions 
at the reference time, Vref is the velocity at the reference time, and a^e/ is the acceleration at 
the reference time. Uncertainties in the fit parameters are determined from the covariance 
matrix. Figures [2] and [3] show the polynomial fits for two example stars and the resulting 
values for the kinematic variables for all stars are reported in Table [2l Since the stars' 
motions are assumed to be dominated by the central force from the black hole, we convert 
'^x,ref and CLy^ref into radial and tangential accelerations B All tangential accelerations and 
positive radial accelerations are non-physical and therefore provide a check on the systematic 
errors of the acceleration measurements. Figure H] shows a histogram of the significance of the 
acceleration measurements both in the radial and tangential directions for the young stars 
in our primary sample. While the tangential and positive radial distributions are slightly 
offset (O.Gcr) from zero and broader (l.Scr vs. la) than is expected for a normal distribution, 
any systematic errors appear to impact the results at the < la level. 

The resulting velocity measurements for the young star sample outside the central 
arcsecond are improved by at le ast a factor of 7 when compared with our previous work 



( IGhez et al.lll998l: iLu et al.l 120051) and o ther recently reported Galactic Center proper mo- 



tions (e.g. iGenzel et al.ll2000l : IOttll2003l ). The absolute uncertainties in our proper motions 
are typically ~0.06 mas yr~^ (~2 km s~^), although stars detected in fewer epochs have 
somewhat higher values (0.1 - 0.5 mas yr~^; 4-20 km s~^). Figures [2] and [3] show exam- 
ples of the measurements for two stars in our sample (SO-15 and IRS 16NW), and their 
corresponding proper motion fits with la errorbars. 

In the young star sample, significant (>3cr) acceleration, or curvature, in the plane of 



^This assumption may not hold for stars in a gravitationally bound cluster, such as may be the case for 
the 4 stars in the extended sample that make up the IRS 13 co-moving group; however, the deviations from 
the potential assumed above should result in only 5 — 10% changes in the velocity vectors. 
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the sky is detected only for SO-15 (Figure [2]). This star has the second smallest projected 
separation from Sgr A* in our sample, at p = I'.'O (0.04 pc), and has a projected radial 
acceleration of -0.21 ± 0.05 mas yr~^ or, equivalently, -9.6 ± 2.0 km s~^ yr~^ (see Figure 
E]). SO-15 is more than twice as far from Sgr A*, in projection, than the seven stars with 
previously detected accelerations, which were all within a project e d radius of less than ?4 
(0.016 pc, IChez et al.l[200ol : lEckart et al.ll2002l : fchez et al.l[2005bl : lEisenhauer et al.l[2005[ ). 



The detection of acceleration is important in that it allows us to solve for the line-of-sight 
distance, and thus the three-dimensional position of a star relative to the black hole. For 
a star in the gravitational potential well of a supermassive black hole, the plane-of-the-sky 
acceleration, at a three-dimensional distance r, in cylindrical coordinates is 

_ -GMp _ -GMp 

where p is the plane-of-the-sky radial coordinate and z is the coordinate along the line of 
sight relative to Sgr A*. The magnitude of the line-of-sight distance from Sgr A*, z, can be 
solved f or by adopting a b lack hole mass of M, = 4.4 x 10^ M© and a distance of Ro = 8.0 kpc 



[see ^ iGhez et al.ll2008l ): it is important to note that there is a remaining sign ambiguity 
for z. The resulting line-of-sight distance from Sgr A* for SO-15 is |0.045 ±0.004| pc bringing 
the total separation between SO-15 and Sgr A* to 0.060 pc. 

The remaining stars in our sample have acceleration measurements that constrain the 
line-of-sight distance. While the lower limits of these acceleration magnitudes are not signif- 
icantly different from zero at the 3a level, their upper limits are smaller than the maximum 
allowed acceleration. The maximum possible magnitude of the acceleration for a star at a 
given p occurs when z = 0. When the measured acceleration limits are below this value, 
they provide a lower limit on the star's line-of-sight distance to the SMBH. Figure [6] com- 
pares the measured acceleration limits with the maximum possible acceleration for each star. 
Any 3(7 acceleration limits below the maximum allowed value gives useful constraints on the 
line-of-sight distances. In addition to our explicit measurement for SO-15, our high precision 
astrometric measurements are now yielding 3cr acceleration limits with a median value of 
-0.19 mas yr'^ (-7.3 km s~^ yr^^) that can significantly constrain the line-of-sight distance 
for nine stars in our sample that are located as far as I'.'T (0.07 pc), in projection, away from 
the black hole. 



4. Orbit Analysis 



For a known point-mass Newtonian gravitational potential, a star's orbital elements 
can be fully determined from the measurement of only six kinematic variables. For this 
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analysis, we assume that the central point mass is a black hole with characteristics deter- 
mined by analysis o f the orbit of the star SO-2, which has been observed for nearly one 
complete revolution (lEisenhauer et al.ll2005l : iGhez et al.ll2008l ). Our proper motion analysis 
yields information on five kinematic variables, including two positions, two velocities. 



and one acceleration. The sixth kinematic variable comes from radial velocities measured by 
Paumard et al.l (120061 ) . The reported radial velocities are averaged over several years of obser- 
vations; however, we adopt the same reference epoch, t^ef, as for the proper motion analysis 
since any change in the radial velocity due to acceleration along the line-of-sight should be 



well within the large measurement uncertainties in radial velocity {a^ 



2, re/ 



-20-100 km s 



As described in §3.31 the plane-of-the-sky acceleration can be converted into a line-of-sight 
distance that, when combined with the projected distance, gives the full three-dimensional 
position for a star. Although most of the stars in our sample have plane-of-the-sky acceler- 
ations that are consistent with zero, the upper limits on the magnitude of the acceleration 
provide valuable information by ruling out small line-of-sight distances. We therefore use 
our best-fit accelerations and uncertainties as formal measures of the acceleration when con- 
verting to a line-of-sight distance. Therefore the six measured quantities can be expressed as 
a three-dimensional position and three-dimensional velocity at a certain epoch (tref)- Given 
the properties of the black hole, these kinematic quantities can be translated directly into 6 
standard orbital elements (see Appendix [B|) . 

A Monte Carlo simulation is carried out to transform each star's six measured kinematic 
variables {xref, Href, v^^ref, '^y,ref, v^^ref, o,p,ref) i^to six Orbital parameters (i, fi, a;, e, P, 
To) and their uncertainties. A total of 10^ Monte Carlo trials are run and, in each trial, 
4 + (6 ■ 32) variables are randomly generated; four for the potential parameters and six for 
each of the 32 stars' measured kinematic variables. The four potential parameters are pulled 
from a four- dimensional probability densit y function, PDF(M,, Rq, Xq, yo), based on the 
orbit of SO-2 derived by iGhez et al.l (120081 ). where the black hole's mass and line-of-sight 



distance are centered on M, = 4.4 x 10^ M0and Ro = 8.0 kpc 0, the dynamical center is 
adopted as the origin with Xo and jo defined as zero, and the projected one-dimensional 
probability distributions' RMS errors are [1.0, 1.6] mas for [xo, Uo], 0.3 x 10^ M0 for M,, 
and 0.3 kpc for |^. For each trial, all the stars' orbits are calculated using the same 



^ These values cor respond to a 12-parameter orbit model for SO-2 (i.e. Vz = case) from a n early version 
of Ghez et al. ( 2008h . In this version, local distortions were not corrected ( Ghez et al. 20081 Appendix B); 
but the resulting black hole mass and distance differ by < la from the final reported values. 

^Simulations were also performed using the lower black hole mass and distance reported by 
( Eisenhauer et al. 20051 ). Our results on the detection of only one stellar disk and on the properties of 



the disk are all consistent within la error bars. 
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potential parameters in order to preserve correlations between the potential parameters 
and the orbital parameters such as eccentricity. The kinematic variables for each star are 
sampled from independent gaussian distributions, each of which is centered at the best-fit 
value from Table [2] and has a la width set to the measurement uncertainty. Any correlations 
between the measured kinematic variables are negligible given the small uncertainties in the 
stars' relative angular positions (^0.2%) and velocities (^3%) in the plane-of-the-sky as 
compared to the uncertainties in the black hole mass (~10%) and the accelerations (~60%). 
The distribution for the acceleration, a^, is truncated such that only accelerations of bound 
orbits are alloweci^, which follows from requiring a negative specific orbital energy, 

^ GM 

E = - < 6 

2 r 

and substituting from Eq. [5] to give the acceleration constraint 

'''''' ^ 8(GM)2- 

For each trial and each star, the orbital parameters are computed and the results of all the 
trials are combined into a six- dimensional probability density function (PDF) by dividing up 
parameter space into bins, summing the number of trials in each bin, and then normalizing 
by the total number of trials. This Monte Carlo method is a straight-forward way to combine 
a star's six measurement PDFs and the four- dimensional PDF for the central point mass, 
which shows strong correlations between M, and Ro, to produce a six-dimensional PDF for 
each star's orbital elements, PDF(z, Q, u, e, P, To), which has strong correlations between 
the orbital parameters. The results of these simulations are plotted for an example star, IRS 
16SW, in Figure [7] to show that i and Q are generally well determined and that e, in some 
cases, can be usefully constrained. Similar figures of the orbital parameters for every star are 
shown in Figure Set 7, which is available online in the electronic edition of this manuscript. 

The resulting stellar orbital parameters are constrained by several different factors. 
First, a measured acceleration that is significantly different from zero, such as for SO-15, 
yields the best determined orbit since the line-of-sight distance is confined to a small range 
of values (Figure [H top). Secondly, each star has a maximum allowed acceleration, ap^max = 
I — GM/p^l, at the closest possible distance set by the observed projected radius. Stars 
with measured accelerations more than 3a below the maximum allowed acceleration, such 



^The assumption that the orbits are bound does not effect the results presented in this paper discussed in 
since all unbound orbital solutions yield high inclination (edge-on) orbits and large eccentricities (e l 1). 

Considering only bound orbits simplifies the orbit analysis as we need only consider equations for elliptical 

orbits rather than hyperbolic or parabolic orbits. 
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as IRS 16NW, have strong lower limits on their line-of-sight distances, which translate into 
significant constraints on the direction of the angular momentum vector, L, and can be 
equivalently expressed as constraints on inclination, i, and on the angle to the ascending 
node, Q (Figure [HI middle). Finally, even stars without significant limits on their line-of- 
sight distance from accelerations have some well constrained orbital elements. In particular, 
i and Q are well constrained as a result of the precise measurements for the stellar velocities 
and potential parameters. Furthermore, if the star's total velocity is higher than the circular 
velocity at the two-dimensional projected radius, then it is higher than the circular velocity 
at all distances and only non-zero eccentricity orbits are allowed (Figure [8], bottom). 

The Monte Carlo analysis described above assumes that, in the absence of an accelera- 
tion measurement, the acceleration should be drawn from a uniform probability distribution; 
or, in other words, we adopt a uniform acceleration prior. For those stars that are only in 
the extended sample, the Monte Carlo orbit analysis samples from this uniform acceleration 
prior ranging from the largest allowed acceleration by the projected radius to the smallest 
allowed for the orbit to remain bound. For these stars and for stars in the primary sample 
with acceleration limits that are not significantly smaller than the maximum physically al- 
lowed acceleration, the uniform-a^ prior is an important assumption. To test how sensitive 
our results are to this assumption, we performed the same Monte Carlo analysis as detailed 
above using an alternative assumption that the prior acceleration distribution is uniform 
in z, which shifts the line-of-sight distance PDF to larger values when compared with a 
uniform-Op prior. On a star-by-star basis, the resulting orbital parameters are consistent 
within la for both priors, with one exception. The young star SO-14 has an eccentricity that 
is constrained to be higher than 0.93 {3a) with a uniform-Op prior, while with a uniform-z 
prior, all eccentricities are allowed within 3a. SO-14 is distinguishable from all other stars in 
our sample in that it has a total velocity of only 50 km s~^, as compared to 160-640 km s~^ 
for the rest of the sample. Such a small velocity translates into a very large range of allowed 
line-of-sight distances which are not well sampled by a uniform-Op prior. SO-14's range of i 
and Q are not largely affected by the choice of prior; therefore, we exclude SO-14 from our 
eccentricity analysis, but we keep it in all other orbital analyses. 

To distinguish between these two possible priors, we examine the resulting distribution 
of orbital phases. For a set of stars whose motion is dominated by the supermassive black 
hole and that have been orbiting for more than a few orbital time scales, the distribution 
of orbital phases should be uniform. The distribution of orbital phases for our sample is 
constructed by summing the orbital phase PDFs for all the stars. Figure M shows that while 
the uniform-Op prior produces a population that is uniformly distributed in orbital phase, 
the uniform-z prior produces a distribution that is strongly peaked at (periapse) due to 
the higher occurence of large line-of-sight distances that, for a given velocity, creates an 
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artificial bias towards periapse. Such a strong bias towards periapse is unlikely to occur even 
if some of the young stars reside in a gravitationally bound cluster, such as IRS 13, where 
all the cluster members would have a similar orbital phase. Based on our assumption that 
the distribution of orbital phases should be roughly uniform, we adopt a uniform-a^ prior 
instead of the uniform-z prior in the following sections. 



5. Orbit Results 

5.1. Detection of the Clockwise Disk 

A large number of stars appear to share a common orbital plane based on our analysis, 
which has no prior assumption about the existence of a disk. The orientation of a star's 
orbital plane can be described by a unit vector originating at Sgr A*'s position and pointing 
normal to the orbital plane (n); and, this normal vector's direction can be expressed by the 
inclination angle (i) and the angle to the ascending node using 

sin i cos Q \ 

— sin i sin i7 . (8) 
— cos i j 

The direction of each star's orbital plane normal vector is determined from the joint two- 
dimensional probability density function of i and fi, PDF(i, which is constructed by 
binning the resulting i and f2 values from the Monte Carlo simulation i n a two-dimensiona l 
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Kn, ) 





histogram with equal solid angle bins using the HEALpix framework (jCorski et al.ll2005l ). 
Figure [TU] shows PDF(i, Vt) projected onto the sky as viewed from Sgr A* for the same 
three example stars shown in Figure [HI Figure [TT] shows, for all stars, the contours for the 
68% confidence region, which, on average, covers a sohd angle of SA^ ~0.2 steradian (sr) 
for the primary sample and 0.6 sr for stars found only in the extended sample, which have 
larger proper motion uncertainties. Table [3] & H] list this solid angle, SA^, for each star in 
the primary and extended samples. The bound orbit assumption does not greatly impact 
the size of the SAfi because the orbital parameters i and VL asymptote at large line-of-sight 
distances as can be seen in Figure [H Stars with acceleration limits significantly smaller than 
0'P,max have two isolated solutions because small line-of-sight distances (z) are not permitted 
and at large line-of-sight distances the positive-z and negative-z solutions asymptote to two 
different values of Q (see Figure [H]). Despite this degeneracy, the clockwise (i=90°-180°) 
stars' normal vectors appear to cluster around a common point indicating that many of 
these stars lie on a common orbital plane. 



The directions of the stars' normal vectors show a statistically significant clustering as 



-17- 



measured by the the density of normal vectors in the sky as viewed from Sgr A*. To quantify 
the density of normal vector directions, we use a n earest neighbo r density estimate, which 



is commonly used to identify galaxy clusters (e.g., lDresslerlll980l ). and take the density at 
each point on the sky to be 

_ k 1 

-stars sr (9) 



27r(l — cos^t 



where 6k is the angle to the k^^ nearest star and k is taken to be 6. We calculate the 
expectation value for the density of normal vectors at each point on the sky using the Monte 
Carlo simulation discussed earlier. For each Monte Carlo trial, the sky is divided into 12288 
equal area pixels (0.001 sr) using a HEALpix grid and the density of normal vectors is 
calculated for each pixel. These estimates are then averaged together over all the trials to 
provide an average density per pixel on the sky. The resulting average density of normal 
vectors is nearly the same for a choice of 4th, 5th, or 7th nearest neighbor. Additionally, 
a similar analysis using a fixed aperture to calculate the density of normal vectors at each 
point on the sky produced similar, but less smooth, results as the nearest neighbor approach 
we adopt here. A peak in the density of normal vectors is detected at i = 115° ± 3° and 
Q = 100° ± 3°, which provides direct evidence of a common orbital plane without any prior 
assumptions (see Figure [T2l) . The uncertainty on the peak position is taken as the half- 
width at half-maximum of the peak divided by the square-root of the number of stars that 
are candidate disk members, \/N disk-stars (see below). We also note that an analysis of the 
entire extended sample produces a peak at the exact same position. The mean density of 
normal vectors at the peak is 0.016 stars deg~^ with a negligible uncertainty on the mean 
value (< 10""^ stars deg^^). The significance of the peak is determined by comparing the 
background density of normal vectors, which is defined by the average (0.001 stars deg~^) 
and standard deviation (0.0008 stars deg~^) of all other pixels on the sky after first rejecting 
those pixels (~0.25 sr) that are high outliers (more than three standard deviations). The 
density peak is ~19ct above the observed background density. A second comparison can be 
made to the density expected if the 32 stars in our sample were isotropically distributed 
over 47r steradians. The observed peak in the density is >20 times higher than this isotropic 
density. Thus we conclude that there is a statistically significant common orbital plane of 
young stars. 

The majority of the young stars that are orbiting in the clockwise direction are likely to 
be orbiting in this common plane. A comparison of each star's normal vector to the common 
plane's normal vector allows us to determine which stars are not on the common plane with 
high statistical significance. All other stars are then considered candidate members. First, a 
preliminary estimate of the thickness of the common plane is determined by defining the solid 
angle extent of the plane, SApiane^ encompassed by the contour at which the density drops 
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to half of the peak value. This corresponds to a region with a solid angle of SApiane ~0.1 
sr, which gives a half-opening angle of 0.2 radians (10°) for a cone with the same SApiane- 
Then, each star's probability density function, PDF(z, Q), is integrated over this region to 
determine the probability that the star is a disk member. The orientation of the stars' normal 
vectors have a wide range of uncertainties as expressed by the total solid angle covered by 
each star, so it is necessary to distinguish between those stars that have a low probability due 
to a large n-uncertainty (i.e. large solid angle) vs. those stars that have a low probability 
because they are significantly offset from the common plane. Therefore, we normalize the 
above integrated probability by the probability at the peak of the star's PDF integrated over 
a region that has the same total area as the common plane 



L(not on plane) 
dSA 

plane 



/plane PDF(z,fi) dSA 



/peakPDF(z,(]) dSA 



/peak 

dSA 



(10) 
(11) 



peak 



where SA is the solid angle and L(not on plane) is the likelihood that the star is not on 
the common plane. Those stars with likelihoods, L(not on plane), of greater than 0.9973 
(equivalent to 3a for a gaussian distribution) are flagged as non-members of the common 
plane. The remaining set of stars are considered candidate members of the common plane. 
Table [3] &|4] list [1 - L(not on disk)] for each star and Figure [5] shows the positions of candidate 
members of the common plane in red and non-members in blue. Of the primary sample of 
32 stars, 26 of which are orbiting in a clockwise sense on the plane-of-the-sky, we find that 
22 are possible members of the common plane {N disk-stars = 22). 

The clockwise common plane that we measure is slightly offset from the clockwise planes 
proposed in earlier works. Over-plotted in black on Figure fT2] is the candidate orbita l plane 
proposed by iLevin fc Beloborodovl (120031 ) with updated values from iPaumard et al.l (120061 ) 
for the candidate plane normal vector (solid black) and thickness (dashed black). The 
previously proposed plane was derived by minimizing a statistical metric, K, in order to find 
the best-fit common orbital plane from the velocity vectors of a sample of stars (see Appendix 
lUj) . However, some stars are not members of the common plane and including them in the 
fit biases the result since they have extremely well n aeasured velocities (SO - 15, IR S 16C, 
S3-19). For example, using the K metric approach of iLevin &: Beloborodovl (l2003l ). fitting 
all 26 clockwise stars in our p rimary sample gives i = 128° and Vl = 102° with K = 0.7, which 
is closer to the disk found by lPaumard et al.l (120061 ) at i = 127° and Q = 99°. While fitting 
only the 22 stars that are consistent with the clockwise disk based on our orbit analysis 
gives i = 117° and Q = 98° with K = 0.2. Therefore, using the K metric to determine 
the common plane can produce biased results due to the inclusion of non-members. By 



- 19 - 



combining position, velocity, and acceleration information in order to determine the orbital 
plane for each star, the direction of a common orbital plane can be estimated more robustly. 

The detected common orbital plane is composed of stars dispersed in a disk rather than 
in a single cluster as can be seen from the stars' positions within the common plane shown 
in Figure [131 In this figure, the stars' positions have been converted into a disk coordinate 
system defined as [p, q, n] where n is perpendicular to the disk plane, p is along the line of 
ascending nodes (where the plane of the sky intersects the disk plane), and q = -h x p. For 
each star, all orbital solutions that fall within 10° of the common orbital plane are combined 
to create a probability distribution for the star's position in the disk, PDF(p, g), which is 
shown in Figure [13] {left). Each stars probability distribution is elongated in the g-direction 
due to the large range of line-of-sight distances, z, that are possible within the small range 
of possible disk inclinations for this nearly edge-on plane of the disk. The thickness in the 
p direction is largely set by the uncertainties in the potential parameters (M,, Rq) and 
velocities. The distribution of young stars within the plane shows a range of position angles 
on the plane, consistent with a stellar disk rather than a stellar cluster. 

The CW stellar disk is detected both in our analysis of the primary sample and in a simi- 
lar analysis of the entire extended sample. The additional young stars in the extended sample 
have larger velocity uncertainties and no acceleration information, therefore the Monte Carlo 
orbit analysis samples from a prior probability distribution that is uniform in acceleration 
ranging from the largest allowed by the projected radius to the smallest allowed for the 
orbit to remain bound. We note that even if we ignore the acceleration measurements for 
our primary sample analysis, the CW stellar disk is still detected, although the significance 
is lowered from ~19cr to ~8cr above the background density of normal vectors. Thus the 
additional stars' orbits in the extended sample are still constrained (see Figure [H]), even 
though they have larger uncertainties as compared to the stars in just the primary sample. 
The density of normal vectors from the extended sample analysis shows a peak within 1° of 
the disk's position from the primary sample. 

The analysis of the extended sample shows that ~50% of the young stars reside on the 
CW disk and there is no statistically significant change (> 3a) in the fractional number 
of disk stars at different radii. For reference, the 73 young stars in the extended sample 
are distributed on the plane of the sky with a surface density that decreases with radius 
as p-2.i±o.4 Within a projected radius of 3", the fraction of candidate disk members is 
72% ± 9% (18 out of 25) and at projected radii larger than 3", the fraction of candidate 
disk members is 42% ± 7% (20 out of 48). Given the small number of known young stars, 
Poisson statistics indicate that this change in the fraction of candidate disk members is only 
marginally statistically significant at the 2.6a level. Likewise, the projected surface density 
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for the on-disk and off-disk populations shown no significant difference from each other or 
from that of the total population. Thus the number of candidate disk members does not 
change with radius and roughly half of the young stars reside on the CW disk. 

The K-band luminosity function (KLF) of the young stars does not change significantly 
with radius or when considering stars on and off the disk. To compare the KLF as a function 
of radius, the entire extended sample of young stars is divided into a near sample (r < 3'.'5) 
and a far sample (r > 3'.'5) and the KLF is constructed for each. A two-sample KS test yields 
a probability of 46% that the near and far samples have the same KLF. Similarly, the KLF 
is constructed for stars on and off the disk and a two-sample KS test yields a probability of 
74% that the on-disk and off-disk samples have the same KLF. Finding more young stars 
will allow for a more detailed comparison of the KLF for different subsets within the young 
stars population. 



5.2. Limits on Additional Stellar Disks 



In our primary sample, no common orbital plane is detected for the counter-clockwise 
population of stars; however, our sample is limited to six counter-clockwise orbi ting stars, 
only two of which (IRS 16NE, IRS 16NW) are claimed by iPaumard et al.l (120061 ) to reside 
on the counter-clockwise disk. Out of the 6 counter-clockwise stars in our primary sample, 
we find that only IRS 16NE and S2-66 could be consistent with the previously proposed 
counter-clockwise disk. The proposed counter-clockwise disk may have a larger radial extent 
than is covered by our observations, so in order to fully explore whether our lack of detection 
of a 2nd disk is due to our limited field-of-view, it is necessary to analyze the extended 
sample. As discussed in ^ the uniform acceleration prior adopted for this analysis tends to 
overemphasize f ace-on orbital planes, making it easier to detect the proposed CCW disk, as 
Paumard et al.l (120061 ) suggest it has an inclination of 24°. 



Using the extended sample, our analysis of the density of normal vectors, in the region 
of the proposed counter-clockwise disk, reveals no significant over-density. Of the 73 stars 
in the extended sample, at least 34 are not on the clockwise disk and thus we compare the 
density observed in the region of the proposed counter-clockwise disk to that expected for 
an isotropic distribution of 34 stars. The observed density of normal vectors in the region of 
the counter-clockwise disk is 2.4x10"^ stars deg~^, which is only a factor of 3 above what is 
expected for an isotropic distribution and is less than la above the background over the rest 
of the sky (excluding the clockwise peak). This density of normal vectors corresponds to only 
3 stars within 1 9° of the putative CCW disk, where 19° is the disk thickness proposed by 
Paumard et al.l (120061 ) . and is consistent with random fluctuations of an isotropic distribution 
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having the n-uncertainties shown in Figure [TH We estimate that this analysis is capable 
of revealing, at the 3cr level, a stellar disk with more than 7 stars within a solid angle cone 
of radius = 19° at the location of t he proposed CCW di sk; thus the proposed CCW disk 



containing 17 stars as suggested by iPaumard et al.l (120061 ) should have been detected with 
this approach. 

There are several principle differences between our analysis and that in earlier works. 
First, previous works make the a priori assumption that a disk exists through the use of 
the statistical metric, K, and the results were not compared to a null hypothesis (i.e. no 
disk) to establish the statistical significance of a disk detection. Furthermore, the K metric 
used in previous works suffers from a bias which is described in Appendix O The primary 
goal of our methodology is to minimize the number of a priori assumptions and to fully 
quantify the significance of any disk detected as compared to the null hypothesis that there 
is no disk. Therefore, we choose to search for disks using all the young stars rather than 
first trimming out stars based on projected angular momentum criteria or radii. Also, we 
determine the range of allowed orbital orientations for each star individually rather than 
searching for a disk from a statistical sample of young stars. In this fashion, we utilize not 
only the direction information for a velocity vector, as has been used previously, but also the 
physical relationships between the magnitude of the velocity and the positional information. 
This method allows for no disk to be detected, while the previously used statistical tests 
assumed a disk model and, therefore, must be compared to the no-disk hypothesis using 
simulations of isotropic populations. Without the simulations, the significance of any disk 
detection via the K metric cannot be fully quantified. The resulting distribution of orbits 
from our analysis is consistent with the hypothesis of a single, clockwise disk plus a more 
randomly distributed population. 



5.3. Properties of the Clockwise Disk 

We now examine, in detail, the properties of the detected clockwise disk. With the 
identification of a single stellar disk and a candidate list of disk members, we investigate the 
following: (1) the thickness of the disk, (2) the radial profile of the disk, (3) the azimuthal 
isotropy of the disk, (4) the eccentricities of stars in the disk, and (5) the luminosity function 
of the stars in the disk. These properties are critical for distinguishing between in situ and 
infalling cluster formation scenarios, as well as for understanding the dynamical evolution of 
the young stars both on and off the disk. 

The observed disk of young stars has a significant intrinsic thickness; however, the 
vertical velocity dispersion is less than previously determined. To measure the thickness of 



-22- 



the disk, the dispersion of the velocities out of the plane (along the n direction) is calculated 
from all candidate disk members by projecting each star's three-dimensional velocity vector 
along the disk's normal vector to give f^. The measurement uncertainties in both v and n 
are propagated through this coordinate transformation. The intrinsic velocity dispersion is 
calculated using 



cr, 



n,intrinsic 
n,mtrinsic 



^n,measured ^n,bias 



disk— start 




(12) 

^ disk — star s disk — star s \ 

XI - error^{vri,i) (13) 



i=0 



i=0 



where the bias term, ais,bias, is 19 km s"^ and accounts for added dispersion as a result of 
uncertainties in the measurements. The resulting intrinsic velocity dispersion is 28 ± 6 km 
s~^, which is significantly different from zero, thus a finite thickness is required. However, this 
velocity dispersi on is a factor o f 2 srn aller than that found using the previously proposed 
disk solution of iPaumard et al.l (120061 ) and is slightly smaller than the value reported in 
Beloborodov et al.l ([20061) due to our improved identification of candidate disk members. The 
disk's thickness can be expressed as the ratio of the vertical scale heig ht to radius, h/r 
o'uMtrinsic/ < \v\ >, and is 0.08 ± 0.02. Following a similar analysis to iBeloborodov et al. 



( 120061 ) ■ but with the above relationship between h/r and the velocity dispersion, the disk 
thickness can also be described using a gaussian distribution of inclination angles about the 
disk plane with a standard deviation of A6 and is related to the scale height of the disk by 
h/r ^J\/2[^Q. This yields a dispersion angle of = 7°±2° for the young stellar disk. This 
more rigorous determination of the disk thickness is consistent with the thickness we derived 
in §5. II from the half-width at half-maximum of the peak in the density of normal vectors; thus 
the selection of the candidate disk members is likely robust. In comparison, the previously 
proposed dis k solutions yield a di sk th ickness of h/r = 0.2 (A ^ = 14°) and h/r = 0.1 
(A^ = 9°) for lPaumard et al.l (120061 ) and IBeloborodov et al.l (120061 ) . respectively. We caution 
that all of these conversions from velocity dispersion to disk scale height and dispersion angle 
assume circular orbits and an isothermal disk structure. From our analysis, we note that the 
out-of-the-plane velocity dispersion shows no statistically significant variation with radius in 
the disk both for the primary (difference is la ~ 7 km s~^) and the full extended samples 
(difference is la ~ 14 km s"^). Therefore, the observations are consistent with a thin disk 
of uniform velocity dispersion at all radii. 

The surface density of stars in the disk falls off rapidly as a function of radius. In order 
to extend the radial coverage, we consider the entire extended sample in this analysis. The 
young stars that are candidate disk members have constraints on their three-dimensional 
radii if we limit their orbital solutions to those close to the disk plane. Thus the disk's 
surface density can be determined as a function of three-dimensional radius rather than just 
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the projected two-dimensional radius as discussed at the end of §5.11 The distribution for 
each star's position within the disk plane, PDF(p, g), is constructed from orbits that are 
within 10° of the disk and is shown in Figure [T3i Then the disk's surface density at each 
radius is computed numerically by sampling the PDF(p, q) 10^ times for all the candidate 
disk members and constructing a radial histogram for each trial. The radial histograms are 
combined for all the trials to find the peak and 68% confidence bounds for the expected 
number of stars at each radius. This is converted into an azimuthally integrated surface 
density by dividing by the area of a ring at each radius. This method of constructing the 
surface density captures both the measurement error in the individual stars and the finite 
thickness of the disk, which has not been incorporated into previous estimates. The resulting 
azimuthally averaged surface density on the disk is shown for the extended sample in Figure 
[T^and has a best-fit p ower-law profile of r~2 '^='=°-^. This is consistent with the previous results 
(jPaumard et al. 2006 ). but our analysis accounts for the uncertainty in each stars line-of- 



sight distance due to the finite disk thickness and, therefore yields a larger uncertainty on 
the power-law index. 

Visual examination of the stars' positions in the disk plane (Figure [T3i) suggests there 
may be some anisotropy as evidenced by the clustering of stars on the lower part of the disk; 
however, this over-density is only marginally statistically significant based on the following 
analysis. In order to search for non-uniformities, we compare the observed stellar surface 
density of the extended sample within the disk plane with the surface density expected for an 
azimuthally symmetric disk. The observed stellar surface density is measured by sampling 
from all stars' PDF(p, q) for 10^ trials and calculating the stellar surface density over a grid of 
points in the disk plane for each trial. For each point on the disk plane, the surface densities 
from all trials are combined, yielding the most probable surface density with uncertainties. 
The resulting two-dimensional map of observed surface densities is then compared to the 
expected surface densities for an azimuthally symmetric disk by subtracting the two values 
and dividing by the uncertainties. This produces a surface density excess map that shows 
the significance of any excess. The disk shows a marginally significant (~ 3a) over-density 
on the front side (q < 0) of the disk and a corresponding under-density on the back side (q 
>0). 

A few candidate disk stars show evidence for eccentric orbits. To determine whether any 
of the stars' eccentricities are consistent with a circular orbit, the six-dimensional probability 
density function for the orbital parameters is marginalized and re-expressed as a PDF for 
the eccentricity vector (see Appendix [Bl) . PDF(ea;, e^, e^). The magnitude of this vector 
is the orbital eccentricity and the direction points along the semi-major axis towards the 
periapse position. The PDF for the eccentricity vector cannot be further marginalized to 
produce a PDF of the eccentricity magnitude without introducing a bias due to the positive. 
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definite nature of a vector magnitude. This is the same bias term as described in the velocity 
dispersion analysis; however, unlike the velocities, the eccentricity distributions are strongly 
non-gaussian and the bias term cannot be easily accounted for in the marginalization. The 
peak of PDF^e^, Cy, e^) gives the unbiased orbital eccentricity and the 99.7% confidence 
interval of the three-dimensional distribution is used to determine the range for the one- 
dimensional eccentricity. Tables [3] and |4] show the 99.7% confidence range of the eccentricities 
for all stars in the primary and extended samples. Also, Figure [16] shows the eccentricity 
99.7% confidence lower limit for the candidate disk members in red, non-disk members 
in blue, and excludes SO-14 (see [g]). When considering all possible orbital solutions, the 
resulting eccentricity ranges show that 2 candidate disk members from the primary sample 
have 99.7% confidence eccentricity lower limits of greater than 0.2. Restricting the possible 
orbital solutions to only those having normal vectors oriented within 10° of the disk normal 
vector increases the number to 8 candidate disk members with 99.7% confidence eccentricity 
lower limits larger than 0.2. 



We find high-eccentricity stars in the disk, similar to the analysis of iBeloborodov et al. 



(120061) in which they assumed an infinitely thin disk. However, our analysis incorporates the 
finite thickness of the disk and places statistical errors on the eccentricities for individual 
stars. 

The average eccentricity of the entire population is not yet well constrained. The ec- 
centricity for the stellar disk is determined using the eccentricity vector. For each candidate 
disk member, orbital solutions are selected whose normal vectors point within 10° of the disk 
normal vector. These orbital solutions are combined for all the disk stars by averaging their 
PDFs to create a combined probability distribution for all stars' eccentricity vectors, which 
is then projected into the disk plane and plotted in two dimensions (Figure [T7I). This two- 
dimensional probability distribution gives an unbiased estimate of the eccentricity magnitude 
and shows that while the characteristic disk eccentricity peaks at e=0.22, it is consistent with 
e=0.0 — 0.8 at the la level, reflecting the large eccentricity uncertainties for the majority of 
the candidate disk members. 



6. Discussion 

The kinematic analysis of the young stars in the central parsec around our Galaxy's 
supermassive black hole has implications for the recent star formation history in this region. 
Our first attempt at determining individual orbits for young stars that reside outside the 
ce ntral arcsecond shows d e finitiy e evidence for the clockwise-rotatin g disk that was sugg ested 



by iLevin &: Beloborodovl (120031 ) and was subsequently refined by iGenzel et al.l ( 120031 ) and 



- 25 - 



Paumard et al.l (120061 ) . Our results do not show a statistically significant second disk. The 
presence of a single stellar disk eliminates the need to invoke two distinct starburst events 
occuring roughly 6 Myr ago and greatly simplifies the demands on both in situ and inf ailing 
cluster scenarios. For instance, in the self-gravitating gas disk scenario, the detection of only 
a single stellar disk hfts the requirement for a second disk to rapidly build up gas, fragment, 
and form stars within 1-2 Myr of the formation of the first disk. Likewise, for the infalling 
cluster scenario, the presence of only one stellar disk means that the frequency of such infall 
events is half that required for the existence of two disks. On the strength of our confirming 
only one stellar disk, we consider whether all of the young stars within the central parsec 
may have formed in a single burst of star formation. 

Such a scenario must explain not only the observed clockwise stellar disk, oriented at 
i ~115° and Q ~100°, but also the presence of roughly half of the young stars from our 
extended sample on more isotropically distributed orbits out of the disk. In the single 
starburst scenario, the out-of-the-disk stars could either be generated during the formation 
process or could intially be in the disk and then perturbed through subsequent dynamical 
evolution. Self-rel axation of the disk has not had suffic i ent ti me to produce the out-of-the- 
plane population (lAlexander et al.l 120071 : ICuadra et al.l 120081 ) . b ut other mech anisms have 
been proposed such as scattering by an inward- migrating IMBH (|Yu et al.ll2007l ). Currently, 
our results show that the on-disk and off-disk populations of young stars look very similar 
outside the central arcsecond (0.04 pc) both in terms of the K-band luminosity function and 
the surface density profiles that decreases at larger projected radii as oc r^^. However, the 
number of young stars in the disk drops at radii smalle r than 0.08 pc; and at radii of 5,0-04 



pc, n one of the observed young S-stars are in the disk (iGhez et al.ll2005bl : lEisenhauer et al. 



20051 ) . This drop in the number of disk stars at small radii may be the result of resonant 
relaxation or other dy namical processes if the cent ral arcsecond S-stars are a continuation 
of the disk population (IHopman fc Alexanderll2006l ). Thus, if dynamical evolution produced 
the off-disk population, then the dynamical process must not be a strong function of radius 
beyond 0.08 pc. 

Our distributions show that a potential problem with the single starburst scenario is the 
presence of the apparent massive star cluster, IRS 13, located ~4" from the supermassive 
black hole (jMaillard et al.ll2004l : ISchodel et al.ll2005l ). The cluster's orbit is not in the disk 
plane and, given the proposed mass of IRS 13 (>10^ Mq ), it is unlikely that it could have 
been ejected from the disk. However, the definition of IRS 13 as a cluster and the derived 
mass is based on observations of only 3-4 bright stars and is complicated by enhanced dust 
and gas emission in the vicinity. More data are needed to determine the total mass of IRS 
13 and its relationship to the disk stars. 
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Our results also have implications for the star formation mechanism. For both infalling 
cluster and in situ formation scenarios, we consider whether the observed characteristics of 
the young stellar disk can be explained. We observe a stellar disk with an out-of-the-disk 
velocity dispersion of 28 ± 6 km s~^. Additionally, if we consider only orbital solutions 
within the disk (disk prior), we find that at least 8 of the 22 candidate disk stars have 99.7% 
confidence lower limits on the eccentricity of greater than 0.2. Therefore, any formation 
scenario should explain not only a single thin stellar disk but also allow for non-circular 
stellar orbits of some stars in the disk. 

First, for the infalling star cluster formation scenario, some of the disk properties we 
observe are well explained and others appear difficult to reconcile with this model. For 
instance, eccentric orbits are easily produced. Stars that are stripped from a cluster as it 
spirals in should have a similar inclination and eccentricity as the cluster itself. Therefore, 
an infalling clus ter with an initially ecceri tric orbit will produce a disk of stars with similarly 
eccentric orbits (IBerukoff fc Hansenll2006r). Previous stu dies have observed co -moving clumps 



of stars, such as IRS 16SW (iLu et al 



20051 ) and IRS 13 (ISchodel et al.ll2005f ). that appeared 



to support the infalling cluster formation scenario as they could be the remnant core of 
the dissipated cluster. We tentatively observe evidence for an over-density of stars on the 
front half of the disk at the position of the IRS 16SW co-moving group. However this over- 
density may be explained by the effects of extinction that reduces the number of young stars 
identified on the back half of the disk at a given magnitude. The extinction is highly variable 
throughout the reg ion and the back half of the disk is beh ind a patch of higher exctinction 
(AAx = 0.3 - 1.4; IScoville et al.ll2003l : ISchodel et al.ll2007l ). Thus the apparent overdensity 
on the front half of the disk, corresponding to the IRS 16SW co-moving group, can perhaps 
be ascribed to differential extinction. More data are needed to confirm the observed disk 
asymmetry and to determine whether the cause is extinctio n. Our results yield a steep 
radial profile for the young stars in the disk, as also found by iPaumard et al.l (120061 ) . which 
appears to be inconsiste nt with the fiatter profile expected for an infalling cluster (r~°'^^, 
Berukoff &: Hansenll2006l ). We note that mass segregation is observed in massive star clusters 



that are only a f ew million years old (IHillenbrand fc Hartmanru Il998l : iFischer et al.l Il998 



Stolte et al.ll2006l ). Any mass segregation that existed prior to the cluster's dissolution may 
impact the observed radial profile as the massive stars would have resided preferentially in 
the cluster core and would therefore have been deposited at the smallest radii. Thus, the 
massive stars O stars that we observe today may have a steeper radial profile than the entire 
young star population. Additiona lly, the lack of X-ray emission from pre- main-sequence 
stars (INayakshin fc SunyaevI |2005| ) is not well explained by an infalling cluster model. A 
larger and deeper survey for young stars over the central ~5 pc could definitively rule out 
this scenario if the tidal tails of the disrupted clusters are not detected. 
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Some theories of in situ star formation take place in a circular gas disk. Such a gas disk 
can be built up from a steady inward migration of material or from many small cloud-infall 
events and the disk would circularize prior to becoming massive enough to form stars from 
self-gravity (> 10*^ Mq ). Such a formation scenario would most likely produce a steep radial 
profile in agreement with our observations. Our observations of over 30% of the candidate 
disk members with eccentricities greater than 0.2 appears to be inconsistent with an initially 
circular disk of stars and a normal initial mass function. A disk of stars on initially circular 
orbits and with a normal IMF will relax over 6 My r and produce a therrn al distribution of 
eccentricities with an rms eccentricity of 0.15 or less (lAlexander et al.ll2007l ). For such a disk, 
only 4 out of 22 stars should have eccentricities higher than 0.2, compared with the 8 out of 
22 observed when a disk prior is imposed on the primary sample. Therefore, in order for the 
disk to have been initially circular with a normal-IMF, some additional dynamical processing 
other than self-relaxation is needed. Other possibilities are that the initial mass function 
may have been top-heavy, the binary fraction may have been extremely high, or IMBHs 
could have formed, all resulting in faster relaxation to higher eccen tricities, but these are 



not sufficient to exp lain the out-of-the-disk population of young stars ( Alexander et al 



Cuadra et al. 



Cuadra et al. 



20081) . The gas disk formation scenario may be modified (lAlexander et al. 



2007; 



2007; 



20081 ) to accommodate the observed high stellar eccentricities and out-of-the 



plane population by building up a massive gas disk in a si ngle cloud infall or a cloud- cloud 
collis i on event, in which the clouds are on eccentric orbits (jSanderd Il998l : IVoUmer fc Duschl 



200ll : iNayakshin et al.l 120071 ). The gas disk w ould then have a high eccentricity for a shor t 



period of time during which stars might form (lAlexander et al.ll2008l : iBonnell fc Ricdl2008l ). 
The cloud-cloud collision scenario may yield both a thin stellar disk and a more distributed 
population of stars at larger radii with a range of angular momenta as a result of the complex 
interactions and shocks during the collision. It is also conceivable that a cloud-cloud collision 
scenario might give rise to out-of-the-disk clumps of gas that could form a cluster such as 
IRS 13. Refined estimates of the eccentricity and inclination distributions of the young stars 
and more detailed theoretical analysis are needed to investigate the viability of this scenario. 



7. Conclusions 

In summary, the advent of laser guide star adaptive optics has allowed us to retroactively 
improve our 11 year astrometric data set used for monitoring stars orbiting our Galactic 
Center. This has increased our proper motion precision, with resulting uncertainties of ~3 
km s~^, and allowed us, for the first time, to make measurements of and place limits on 
accelerations for stars outside the central arcsecond out to a radius of 3'.'5, with typical 3a 
acceleration limits of -0.19 mas yr~^. By combining our improved stellar positions and proper 
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motions with radial velocity information from the literature, we compute orbits for individual 
young stars proposed to lie in stellar disks orbiting the supermassive black hole. The orbits 
for the young stars confirm only a single disk of young stars at a high inclination rotating 
in a clockwise sense and there is no statistically significant evidence for a second disk. Stars 
within the well-defined, clockwise disk have an out-of-the-disk velocity dispersion of 28 ± 
6 km s~^ and several stars have high eccentricities. These disk properties suggest that star 
formation may have occurred in a single event, rather than the two events previously needed 
to explain two stellar disks; however, there are open questions as to how ~50% of all young 
stars can be perturbed out of the disk plane and whether the apparent compact cluster, IRS 
13, which is not part of the stellar disk, requires a separate star formation or dynamical event. 
Future directions include (1) obtaining new LGSAO data sets with improved astrometry to 
measure accelerations for the young stars at all radii and (2) identifying new young stars 
within the central parsec in order to better constrain the orbital properties of these stars 
and to study in detail the distribution of eccentricities and semi-major axes for stars both 
in and out of the disk. 
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A. NIRC Speckle Distortion 

In the speckle data sets, optical distortions, introduced by the NIRC reimager, are small 
near the center of the field-of-view where Sgr A* was positioned, but grow to dominate the 
positional uncertainties for stars located more than ~0'.'5 from Sgr A* (see Figure [18] and 
§3.ip . Now, utilizing images of the Galactic Ce nter obtained wit h NIRC2, which has optical 



distortions characterized at the ~2 mas level (IGhez et al.ll2008l ). we can, for the first time 



similarly quantify and correct the optical distortions in the NIRC reimager speckle data sets. 

Images of the Galactic Center were obtained with both NIRC and NIRC2 on consecutive 
nights during July 2004 and the NIRC2 images were used as a reference coordinate system. 
The individual NIRC speckle exposure times are only 0.1 seconds and have insufficient signal- 
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to-noise to detect more than the brightest 5 stars. Exposures were obtained in sets of 100 
and each set is combined to produce a single image in which approximatly 100 stars are 
detected. It is assumed that the images are mostly stationary on the NIRC detector during 
each set of exposures. For each stacked image, the stars' positions are compared to those 
in the NIRC2 image and the offsets are mapped into NIRC detector coordinates (see Figure 
[T9l left). In this fashion, a distortion map is built up from many stacks of images which are 
dithered and rotated such that stars fall on many different positions on the detector. The 
distortion solution was obtained by fitting the distortion map with polynomials of the form 

{x' + 128) = ao + ai{x - 128) + a2{y - 128) (Al) 

{y' + 128) =bo + h{x- 128) + b2{y - 128) (A2) 

where the best-fit distortion parameters are listed in Table O The new distortion solution 
improves the RMS residual errors per stack by a factof of 3 to ~3 mas (Figure 1201) . which is 
further reduced in the final image by averaging the dithered stacks. Higher-order polynomial 
terms did not sufficiently improve the fit to warrant inclusion. The above solution is applied 
after the initial application of the standard NIRC distortion correction. The map of positional 
differences between stars in the NIRC and NIRC2 images before and after the NIRC-reimager 
distortion correction is shown in Figure [T^ (right). The resulting radial dependence on the 
RMS positional uncertainty is greatly improved and is shown in Figure [181 which plots many 
stars' RMS residual offset from their best-fit proper motions across all epochs. In the final 
analysis of the speckle data, the relative astrometric uncertainty is ~2 mas. 



B. Analytic Orbit Equations 

The orbit of a star in a known point source potential can be derived from a single 
measurement of a star's orbital state vector. At epoch t^e/, the orbital state vector is 
usually described by the star's position, r, and velocity, v, relative to the central mass. For 
the analysis in this paper, the state vector is estimated using measurements of the three- 
dimensional velocity, v = [v^, Vy, f^], and the projected position, r2D = [x, y], and z is derived 
from the radial acceleration on the plane of the sky. For brevity, we have removed the ref 
subscript notation and all of the above variables are measured at tref- Orbital trajectories 
are then inferred from conservation of energy, specific angular momentum, and eccentricity 
(e, h, e), which are related hj e-h = and |ep — 1 = 2 e h'^/GM giving 5 constants of motion 
plus an undetermined reference time. Equivalently, the orbital trajectory can be expressed 
using the standard Keplerian orbital elements: period (P), eccentricity (e), time of periapse 
passage (To), inclination (z), position angle of the ascendin g node (Q) , and t he longitude 
of periapse {to; see Equations [B5l EH EH EH EH EH and 



Ghez et al. 2005b. for detailed 
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descriptions of these orbital parameters). The 3D position and velocity state vectors can 
be used to calculate the orbit of the star around the black hole (by algebraic manipulation 
of Kepler's Laws). Here we present the analytic expressions used to compute the orbital 
elements from the state vectors. 

Orbit determination for the young stars in our sample is tractable because the mass and 
position of the black hole are determined by independent means, namely the well determined 
orbits of stars much closer to the black hole. The coordinate system is set such that Sgr A* 
resides at the origin, x and y increase with right ascension and declination, and z increases 
with the line-of-sight distance from the Earth to Sgr A* with z=0 at the location of the black 
hole. Combining the two state vectors, r and v, and the black hole mass, there are three 
intermediate vectors that describe the geometry of the orbit both in three-dimensions and 
projected onto the plane of the sky. These are (1) the specific angular momentum vector, h, 
which points normal to the plane of the orbit, (2) the eccentricity vector, e, which points in 
the direction of periapse, and (3) the ascending node vector, Q, which points to where the 
star passes through the plane of the sky moving away from us, and are given by 

h ^ fxv (Bl) 

_ V X h f , 

' = GM "M '^'^ 
n = hxz. (B3) 

The semi-major axis can also be calculated as an intermediate quantity 

Then the five standard orbital parameters that describe the shape and period of the orbit 
are then 

i — arccos ( \ ] (B5) 

e = \e\ (B6) 

((z X h) 
— f—^ I (if e • z < then u; — 27: — lo) (B7) 
I Z X 

Q. X 

Q = arctan I ^— ) (B8) 



[yr]J V \[AU]J V M 



(B9) 
(BIO) 
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where i = if the orbit is in the plane of the sky and fl is measured East (ct) of North [y). 
The remaining orbital parameter is the epoch of periapse passage and can be computed in a 
number of different ways. We first compute several intermediate quantities of interest such 
as the Thiele-Innes constants (A,B,C,F,G,H)) and the eccentric anomaly as shown below: 

A = a(costi;cosi7 — sintusinficosi) (Bll) 

B = a(cosco'sinf2 + sincucosi^cosi) (B12) 

F = a(— sin cos i7 — coscij sin i7 cos z) (B13) 

G = a(— sincj sin f2 + cosa; cosf2 cosi) (B14) 

„ Arrc — BVy 1 

sinE = — ^ ^ , B16 

AG - BF v^r^^ ^ ' 

^ ( sin£'\ ,^ 

And finally the epoch of periapse passage are calculated from these intermediate quantities 
using 

To = tref-^{E-esmE). (B18) 



C. K Metric 



The previously proposed planes were derived by minimizing a metric that lLevin fc Beloborodov 



(120031 ) call x^, but we call K, and which is defined as 



- (ci) 



K = 



where is the number of stars, Vi is the velocity of each star, the velocity 

uncertainties for each star, and n is the normal vector to the disk plane that is found in the 
fitting process. This metric is used to find, statistically, the best-fit common orbital plane 
from the velocity vectors of a sample of stars. The K metric suffers from several shortcomings. 
First, the K metric is described as a metric; however, standard minimization takes 
the form of (data - model)^/(data errors)^ where the data errors have no dependency on 
the model parameters. The K metric includes the model parameters in the data-error term 
and does not necessarily have an expectation value of 1 for normal errors. The appropriate 
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function to minimize in order to find the best-fit common orbital plane can be derived from 
maximum likelihood theory if we assume that the likelihood function is given by 

^ 1 

where cjj depends on the disk model parameters that are being sought by 

(Ti = + Ka^^J^ + {n.a^^J^. (C3) 

Standard practice is then to take the logarithm of the likelihood, L, and minimize the 
resulting function in Equation IC5I in order to find the best fit disk model parameters. The 
above likelihood function then becomes 

InL = -^ln(2vr)-^lna, + 5^-^^ (C4) 

1=1 i=l * 

-21nL = Arin(27r) + 2^1na, + ^ ^""'^'^ (C5) 

1=1 i=i ^« 

and the first two terms are constant and do not factor into finding an extremum in the above 
equation. The third term on the right-hand side is the K metric previously used to determing 
the disk parameters. However, the second term on the right-hand side also depends on the 
free parameters in n and must be included in the minimization process. This extra term 
that has not previously been included in the disk fitting process has the full form 

In ^ {n^fyv^ -)"^ + K(^,;,,J^ + {n^(^v, ;)'^ (C6) 

and standard chi-squared probability functions cannot be applied. Second, even when ac- 
counting for the extra term, the metric can still introduce substantial bias. In particular, 
radial velocity uncertainties, cr^^ ^ , are larger than the proper motion errors by a factor of 
2 on average in previous publications. During K-minimization, this over-weights solutions 
with a larger rtz resulting in a bias against edge-on planes. Finally, in order to properly 
evaluate the probability of obtaining a given value of the K-metric by random chance, one 
must perform simulations of an isotropic distribution of stars. However, such simulations 
are extremely sensitive to the input distribution of semi-major axes and eccentricities which 
are not yet well constrained by observations. Thus, when utilizing such statistical tests for 
finding a common orbital plane, it is difficult to compare to the null hypothesis - an isotropic 
distribution of stars - and to quantify the significance of a disk. 
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''Dates are computed as the weighted average of UT dates from the individual exposures. 

''Filters used include K (Ao=2.2 ^m, AA=0.4 fim), K' (Ao=2.12 /xm, AA=0.35 /^m), Kco (Ao=2.289 /^m, AA=0.048 fim), and 
Kcont (Ao=2.270 fim, AA=0.030 fim). 

'^The number of frames used in the final combined image. 

''The turnover of the number of stars at a given magnitude provides a rough estimate of the completeness limit. 

°The average positional uncertainty due to centroiding in each epoch is estimated from a set of 25 stars detected in all epochs 
and brighter than K<13. The two LGSAO epochs positional errors include and additional term of 0.88 mas to account for residual 
distortion. 

*jData ori ginally reporte d in (plOhez et al.l lll998l ') (2) \Ghez et "akl l l2000h ('3') lGhez et "all ll2005d) (4) \Lu et al]||2005f) (5') lGhez et"al] 
l l2005al ') fe) iRafelski et al.l ||2007|) . 

^Five exposures were taken in each of two narrow-band filters with different exposure times, but similar sensitivity and astrometric 
precision. All frames from both filters were combined in order to extract astrometric measurements from this data set. 



Table 2. Proper Motions for Young Stars 



Name 


K 


^epochs 


Epoch 


Radius 


ARA ^ 


ADEC ^ 




^dec 


Vz ^ 


dp 




O-tan 


AltName ^ 




(mag) 




(year) 


(arcscc) 


(arcsec) 


(arcsec) 


(mas/yr) 


(mas/yr) 


(km/s) 


(mas/yr^) 


(mas/yr^) 




oU-i4 


1 Q ? 


zy 


onni oon 
zuui.zyu 


n fio 

U.oZ 


n 77n 

-U. ( (U 


n ovn 

-U.Z ( U 


1 fio -U n nfi 

l.DZ it U.UD 


n /ifi -U n 07 

-U.4D it U.U( 


1 /I _|_ ACX 
-14 It 4U 


0.05 i 


0.05 


n no -I- n oc; 
U.uy it U.UO 


T?1 A 

rji4 


OU- 10 


1 7 
lo. < 


OQ 
zy 


onn 1 fisn 

ZUUl.DoU 


n 07 
u.y I 


n OQO 
-u.you 


U.ZoU 


QO -1- n 07 
-o.oz m u.u i 


1 n 9Q -1- n OS 

-lU.ZO It U.Uo 


A94 -i- 70 
-4Z'y! It / U 


-0.21 i 


0.05 


n 1 n _|_ o na 
U.IU It U.UD 


"PI fl 
HjID 


Q 1 Q 


1 O 1 

lz.± 


on 

zy 


onni riQA 
zUUi.yoU 


n QQ 
u.yo 


n A An 
U.44U 


n Q7n 
U.o ( y 


1 Q QQ _L n nc; 
-lo.oo it U.UO 


o nn _L n 
z.UU it U.Uo 


fiQ -1- An 
Do it 4U 


-0.10 it 


0.03 


n nQ J- o HQ 
-U.Uo it U.Uo 


T?i c; 


Ql 
oi-z 


1 /I Q 

14. y 


OK 
ZD 


oom nan 

ZUUl.oDU 


1 m 

l.Ul 


n noc^ 

-U.UZO 


1 nn7 

-l.UU/ 


1 1 7n _i_ n 1 Q 
11. / U it U.lo 


n UK _i_ n 1 Q 
-U.DO it U.lo 


OH -I- QO 

ZO it oU 


0.01 it 


0.08 


n 1 7 -U n riQ 
U.l 1 it U.Uo 


T71 7 


Ql « 


14. Z 


OQ 

zy 


onm AQfi 

ZUUl.OcSU 


1 nc 

l.Uo 


-U.DOl 


-U.OOO 


7 dA -U n 1 n 
/ .04 ni U.iu 


A KQ -1- n 1 n 

-4. Do it U.IU 


-oD4 it 4U 


0.04 it 


0.06 


n na -i- n na 
U.UD it U.UD 


TTI Q 
HflO 


TT? Q 1 filVRA/' 
Irto iOIN W 


lU.l 


OQ 

zy 


onm Kfin 

ZUUi.ODU 


1 OO 

1 .zz 


n noQ 
u.uzy 


1 001 
l.ZZl 


fi Qn -U n n^; 

D.OU it U.UO 


n Q7 -U n na 

U.o/ it U.UD 


A A -i- on 

-44 It ZU 


0.04 it 


0.04 


n nc; _i_ n nA 
-U.UO it U.U4 


T?1 Q 


iixo loo 


Q Q 

y.o 


OO 

Zy 


onm f^7n 

ZUUl.O ( U 


1 OQ 
1 . ZO 


1 1 oi 

1 . IZl 


n /107 
u.4y / 


Q 7/1 -U n nf^ 

-o. 1 4 it U.UO 


7 /10 -U n nc; 
1.4Z it U.UO 


1 9f^ -1- QO 
IZO it oU 


-0.04 it 


0.03 


o 1 1 -I- n OQ 

-U.ll ± U.Uo 


iron 


Ol- iZ 


1 Q Q 


Ofi 

Zo 


onni f^nn 

ZUUi.OUU 


1 Qn 
l.oU 


n fiQ7 
-U.oo ( 


1 nnn 

-l.UUU 


o OQ -U n n7 
y.yo it u.u 1 


1 fiQ -U n 07 

-l.oo it U.U( 


9/1 4- QO 
-Z4 it OU 


in 


u.uu 


o nn -U n nfi 
-U.UU it U.UO 


1T01 


Q 1 1/1 


1 O Q 
IZ.O 


OO 

zy 


onm Qcn 


1 QO 

i.oy 


-l.oOO 


n Qno 

-U.oUZ 


/I m -U n nfi 
4.U1 it U.UD 


70 -U n 07 

-D. ( y it u.u ( 


/1Q/1 -U c;0 
-4o4 it OU 


-0.11 it 


0.05 


n no -i~ n nc; 
-U.UZ it U.UO 


1T09 


irto lOo VV 


inn 


OQ 

zy 


onm ziQn 
zuui.4yu 


1 A'i 
1 .40 


l.UOl 


-u.yoD 


D.oU It U.UO 


90 -I- n nfi 

Z.ZZ ^ U.UD 


Q9n -1- A(\ 
OZU It 4U 


-0.08 it 


0.04 


n HQ -1- n (\A 

U.UO It U.U4 


PTOQ 

ilizo 




io.o 


1 T 


onm 1 on 
zUUi.iyU 


1 RQ 


1 aan 

-i.Doy 


n 1 /II 
U.14i 


Q j_ n no 
o.oz it U.Uy 


Q Q/1 -1- n no 

-o.o4 it u.uy 


QAA -L- r;n 
-o44 it OU 


0.03 it 


0.07 


n no _L n n7 
U.Uz it U.U/ 


TTO/I 




10 7 
IZ. 1 


OO 

zy 


onm onn 


1 7n 

1. ( U 


1 IRQ 1 


n /10Q 


a of^ -U n n7 
D.yo it U.U 1 


1 70 _i_ n nc 
-1. ( U it U.Uo 


09/1 -1- c;0 
-ZZ'i ZiZ OU 


-0.09 it 


06 


n nQ -1- n nfi 

-U.Uo zt U.UO 


1!jZ0 




1 1 ^; 

li.O 


OO 

zy 


onm /ion 




n 71 Q 
U. ( io 


1 p;m 
-i.oyi 


1 1 Q _l_ n CiT 

i.lo it U.U ( 


fi Q7 J- n OQ 

-D.o ( it U.Uo 


onfi -1- QO 
zUD ± oU 


0.02 zt 


0.05 


n HQ -1- n op; 
U.Uo ± U.Uo 


Ij^Ofi 
£!jZO 


Q9 A 
oZ-4 


1 O Q 
IZ.O 


OO 

zy 


onm /iQn 


o nv 

Z.U 1 


1 A f^O 
1 .40Z 


1 /I Tfi 
-1.4(0 


AO _i_ n nfi 
D.Dy it U.Uo 


O _|_ n OQ 

Z.Ol it U.Uo 


OQfi -U 90 
ZoD it ZU 


-0.03 it 


0.06 


o ofi 4- n nfi 
-U.UD it U.UO 


TTOQ 


TRQ ifir~^n 

lHo IDOU 


1 n 


zo 


onnn s/in 

ZUUU.04U 


o n7 
Z.U / 


1 QQQ 

i.yyy 


n c^cn 
U.oOU 


1 Qc -U n nfi 
-l.oo it U.UO 


K AQ -i- f\ Ofi 
0.4o it U.UD 


OA 1 _i_ oc: 

Z41 It zo 


-0.02 it 


0.04 


n no -U n nA 
-U.uy it U.U4 


P97 


QO ^1 


1 O 1 
IZ. 1 


OQ 

zy 


onm OQn 
zuui.zyu 


HQ 

z.uy 


1 KQA 

i.oy4 


1 Q/l 
-l.o40 


K sjn -1- n nc^ 
D.oU it U.UO 


1 dfi _i_ n na 
l.DO It U.UD 


oifi -U on 

ZIO It ZU 


-0.02 it 


0.04 


n no -1- n nA 

-U.UZ ± U.U4 


WQn 1 


QO 7 


1/11 
14. 1 


07 
Z ( 


onno Qc^n 

ZUUZ.oOU 


o no 
z.uy 


n Q7Q 

u.y ( y 


1 Q/1Q 


-D.lo it U.ll 


1 m _i_ n 1 1 
l.Ul it u.ll 


QA -\- ^n 
-y4 it ou 


0.04 it 


0.08 


n n7 -U n oq 

-U.U* it U.Uo 


TTOQ 


TOO OQ"\T 


1 n Q 
lU.o 


OQ 

zy 


onm Am 

ZUU1.4iU 




-i.oyo 


1 /lOQ 


oa -U n no 

O.ZD it U.Uo 


A Ai -\- r\ riQ 
-4.41 it U.Uo 


1 QO -U QO 

-lyu it yu 


-0.02 it 


0.06 


n no -U n na 
-U.uy it U.UD 




llxo lOoVV-ilj 


1 1 n 

11. U 


OQ 

zy 


onm A'irx 

ZUU1.4oU 


O 1 7 
Z.l ( 


1 QA« 
l.o40 


1 1/11 
-1 . 141 


A SQ -U n nfi 

4.00 it U.UD 


o OQ -U n ofi 
z.ys it U.UO 


Qfifi _|_ 70 
vSOD it (U 


-0.03 it 


0.04 


o OQ -U n O/I 

-U.UvJ it U.U4 


1TQ9 

1!jOZ 


llxo vJOiN 


11/1 

11.4 


OO 

zy 


onm fiQn 

ZUUl.DoU 


O 1 O 

z. ly 


n n/1 Q 

-U.U4o 


O 1 fiO 

-z. loy 


"1 vo _i_ n 1 o 

1. / Z it U. iZ 


cc 1 K _i_ n 1 9 
-0.10 it U.lz 


fiQ -U 90 
Oo it ZU 


-0.01 it 


0.08 


n 1 Q -U n OQ 

U.lo it U.Uo 


T?QQ 
n/OO 


S2-17 


10.9 


29 


2001.660 


2.26 


1.271 


-1.871 


7.51 di 0.09 


-0.51 ± 0.09 


100 ± 20 
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Note. — All uncertainties are la relative errors and do not include errors in the plate scale, location of Sgr A*, or position angle. 
"Positions as determined from polynomial fitting have relative errors of ~0.4 mas. 
'^Radial velocities and alternate names obtained from Paumard et al. (2006). 
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Table 3. Orbital Eccentricity and Disk Membership 



Name 




1 - L(iiot on disk) 


Eccentricity' 


(All Soluti 


ions) 


Eccentricity (Disk Solutions) 


Direction 




(s|.('ra< iiaiis) 




i'eak 


.'-SfT i^ange 




i*eak 


~>(7 [\ iUl,^( 








CANDIDATE DISK MEMBERS 








S2-16 


U.47 


7.89e-01 


0.60 


0.00 - 


1.00 


0.21 


0.00 - 1.00 


OW 


irsloSW-ilj 


0.18 


6.76e-01 


0.37 


0.00 - 


1.00 


0.37 


0.12 - 1.00 


OW 


81-14 


0.12 


5.80e-01 


0.33 


0.00 - 


1.00 


0.33 


0.13 - 1.00 


cw 


S2-6 


0.15 


5.76e-01 


0.79 


0.17 - 


1.00 


0.30 


0.17 - 0.60 


cw 


bO-0 


0.20 


5.62e-01 


0.64 


0.00 - 


1.00 


0.53 


0.06 - 1.00 




irsiboW 


0.14 


5.38e-01 


0.78 


0.04 - 


1.00 


0.41 


0.29 - 0.91 


U W 


CI in 


0.10 


4.56e-01 


0.41 


0.00 - 


1.00 


0.33 


0.00 - 0.61 


OW 


On A 


0.13 


4.23e-01 


0.69 


0.10 - 


1.00 


0.32 


0.21 - 0.94 


ow 


ol-o 


0.09 


4.06e-01 


0.62 


0.37 - 


1.00 


0.57 


0.45 - 1.00 


ow 


CI o 


0.10 


3.66e-01 


0.26 


0.00 - 


1.00 


0.19 


0.00 - 0.80 


U W 


cn 1 7 


n ni 


3.63C-01 


0.77 


0.00 - 


1.00 


0.40 


0.00 - 0.56 


W 


83-10 


0.08 


3.21C-01 


0.16 


0.00 - 


1.00 


0.67 


0.24 - 0.81 


CW 


cn '7 


0.54 


3.02e-01 


0.76 


0.00 - 


1.00 


0.55 


0.07 - 0.67 


our 


83-25 


0.50 


2.64e-01 


0.76 


0.00 - 


1.00 


0.61 


0.27 - 1.00 


OW 


82-74 


0.24 


2.01e-01 


0.45 


0.00 - 


1.00 


0.15 


0.00 - 1.00 


CW 


81-21 


0.14 


1.67e-01 


0.92 


0.00 - 


1.00 


0.46 


0.04 - 0.79 


cw 


82-19 


0.22 


1.60C-01 


0.57 


0.00 - 


1.00 


0.18 


0.00 - 0.46 


cw 


irsieCC 


0.17 


1.49C-01 


0.62 


0.34 - 


1.00 


0.54 


0.40 - 0.66 


cw 


irs33E 


0.31 


1.42C-01 


0.49 


0.11 - 


1.00 


0.50 


0.27 - 1.00 


cw 


81-3 


0.06 


1.32e-01 


0.34 


0.00 - 


1.00 


0.09 


0.00 - 1.00 


cw 


81-22 


0.27 


1.02e-01 


0.92 


0.00 - 


1.00 


0.68 


0.34 - 0.83 


cw 


80-14 


0.13 


5.35e-02 












cw 






STARS NOT IN THE DISK 








80-15 


0.12 


1.55e-03 


0.30 


0.00 - 


1.00 






cw 


irsieC 


0.07 


5.28e-04 


0.50 


0.00 - 


1.00 






cw 


irs33N 


0.14 


< l.OOe-05 


0.97 


0.00 - 


1.00 






cw 


81-24 


0.06 


< l.OOe-05 


0.98 


0.13 - 


1.00 






ccw 


83-19 


0.57 


< l.OOc-05 


0.80 


0.00 - 


1.00 






cw 


83-30 


0.03 


< l.OOc-05 


0.99 


0.00 - 


1.00 






ccw 


irsl6NE 


0.14 


< l.OOc-05 


0.19 


0.04 - 


1.00 






ccw 


irsl6NW 


0.08 


< l.OOe-05 


0.70 


0.00 - 


1.00 






ccw 


irs29N 


0.01 


< l.OOe-05 


0.99 


0.00 - 


1.00 






ccw 


82-66 


0.53 


< l.OOe-05 


0.97 


0.08 - 


1.00 






ccw 



*No eccentricity is reported for 80-14 since the uniform-ap prior is not appropriate for this very low velocity star. 
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Table 4. Orbital Eccentricity and Disk Membership for Stars Added to the Extended 

Sample 



Name 


SA^ 


1 - L(not on disk) 


Eccentricity (All Solutions) 


Eccentricity (Disk Solutions) 


Direction 




(steradians) 




Peak 


3(T Range 


Peak 


3<T Range 








CANDIDATE DISK MEMBERS 








pa.uniE57 


0.41 


8.57e-01 


0.29 


0.00 - 1.00 


0.34 


o.uo — o.yi 


CW 


irs34W 


0.21 


5.12e-01 


0.20 


0.00 - 1.00 


0.20 


o no 1 on 
O.UO — l.UU 


cw 


paumE72 


1.50 


4.31e-01 


0.81 


0.00 - 1.00 


U.OD 


n no 1 on 
U.UU — i.UU 


CW 


paumE73 


0.98 


2.41e-01 


0.02 


0.00 - 1.00 


U.yD 


n on 1 no 
U.UU — i.UU 


cw 


irs34NW 


0.39 


2.33e-01 


0.04 


0.00 - 1.00 


U.U / 


n on 1 no 
U.UU — i.UU 


cw 


AFNWNW 


1.52 


2.06e-01 


0.99 


0.00 - 1.00 


0.96 


n on 1 no 
U.UU — i.UU 


cw 


paumE69 


0.51 


1.34e-01 


0.02 


0.00 - 1.00 


U.oU 


n 1 o 1 no 
U.iz — i.UU 


cw 


irsQSW 


0.60 


l.lOe-01 


0.05 


0.00 - 1.00 


U.OO 


n on 1 oo 
U.UU — l.UU 


cw 


paumE54 


0.36 


1.03e-01 


0.08 


0.00 - 1.00 


0.17 


O oo O 

U.UU — U.OO 


cw 


irslE 


0.88 


7.35e-02 


0.77 


0.00 - 1.00 


U.yo 


o 1 oo 
U.OO — l.UU 


cw 


irs9W 


0.57 


6.29e-02 


0.02 


0.00 - 1.00 


U.0 1 


O 1 v 1 oo 
U.ii — i.UU 


cw 


paumE87 


0.65 


5.09e-02 


0.09 


0.00 - 1.00 


n 0/1 
u.y4 


O 1 oo 

u.oz — i.UU 


cw 


irslSSW 


0.35 


2.64e-02 


0.02 


0.00 - 1.00 


n o/i 
U.y4 


O Q 1 1 OO 
U.oi — i.UU 


cw 


AF 


0.23 


1.22e-02 


0.11 


0.00 - 1.00 


u.yy 


O 77 1 OO 
U. M — i.UU 


ccw 


irslW 


0.59 


9.61e-03 


0.18 


0.00 - 1.00 


KJ.o 1 


O 1 OO 

U.OO — l.UU 


cw 


irsTSW 


0.06 


8.14e-03 


0.99 


0.00 - 1.00 


u.yo 


O 1 OO 
U.Dtt — l.UU 


cw 


83-26 


0.22 


3.32e-03 


1.00 


0.00 - 1.00 


u.yz 


O 7/1 1 OO 

U. r 4 — l.UU 


cw 






STARS NOT IN THE DISK 








AFNW 


1.08 


1.31e-03 


0.95 


0.00 - 1.00 






cw 


irsl5NE 


1.07 


8.92e-04 


0.99 


0.00 - 1.00 






ccw 


paumE78 


0.61 


4.06e-04 


0.98 


0.00 - 1.00 






ccw 


irs9SE 


0.89 


7.12C-05 


0.02 


0.00 - 1.00 






ccw 


irs7E2? 


0.56 


6.80e-05 


0.99 


0.00 - 1.00 






cw 


paumE84 


0.22 


4.09e-05 


0.98 


0.00 - 1.00 






cw 


paumESQ 


0.94 


< l.OOe-05 


0.02 


0.00 - 1.00 






ccw 




u.ou 


< l.OOe-05 


0.98 


0.00 - 1.00 






K^\y VV 


paumE86 


0.32 


< l.OOe-05 


0.06 


0.00 - 1.00 






ccw 


paumE82 


0.79 


< l.OOe-05 


0.02 


0.00 - 1.00 






ccw 


paumE75 


0.42 


< l.OOc-05 


0.99 


0.00 - 1.00 






cw 


paumE64 


0.67 


< l.OOc-05 


0.62 


0.00 - 1.00 






cw 


paumE62 


0.58 


< l.OOc-05 


0.02 


0.00 - 1.00 






ccw 


paumE60 


0.21 


< l.OOe-05 


0.99 


0.00 - 1.00 






ccw 


paumE55 


0.61 


< l.OOe-05 


0.95 


0.00 - 1.00 






ccw 


paumE53 


0.49 


< l.OOe-05 


0.97 


0.00 - 1.00 






cw 


paumE52 


0.26 


< l.OOe-05 


0.09 


0.00 - 1.00 






ccw 


paumE42 


0.59 


< l.OOc-05 


0.04 


0.00 - 1.00 






ccw 


irs7W 


0.33 


< l.OOc-05 


0.98 


0.00 - 1.00 






ccw 


irs7SE 


0.81 


< l.OOc-05 


0.07 


0.00 - 1.00 






ccw 


irs7El(ESE) 


1.18 


< l.OOe-05 


0.05 


0.00 - 1.00 






ccw 


irsl3E4 


0.60 


< l.OOe-05 


0.70 


0.00 - 1.00 






ccw 


irsl3E3b 


0.53 


< l.OOc-05 


0.64 


0.00 - 1.00 






ccw 


irsl3E2 


0.59 


< l.OOe-05 


0.60 


0.00 - 1.00 






ccw 
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Table 5. NIRC Reimager Distortion Coefficients 



i 


X (a,) 


Y (b,) 





1.713x10-2 


-2.654x10-2 


1 


9.957x10-1 


-1.759x10-3 


2 


-3.371x10-3 


1.004 
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Fig. Set 7. Orbital Parameters 
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Fig. 1. — Positional uncertainties for stars as a function of stellar brightness (top) and 
distance from the black hole, Sgr A*, which is near the center of the field of view (bottom). 
To show the full range of possible values, the centroiding (solid) and the alignment (dashed) 
uncertainties are shown for the best (1999.559) and worst (1996.485) speckle epochs and one 
of the LGS AO epochs (2004.567). The uncertainties are the median values of all stars within 
magnitude bins of AK = 1 or radius bins of Ar = 0'.'3. Note that alignment uncertainties 
are small compared to centroid uncertainties. 
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Fig. 2. — Measured positions and residuals as a function of time for SO-15, a source with 
a significant, non-zero acceleration measurement, in X and Y (top). Positions are reported 
relative to Sgr A* and do not include the uncertainties in the transformation to the absolute 
coordinate system (i.e. plate scale, position angle, and position of Sgr A*). The best fit 
quadratic polynomial modeling the velocity and acceleration of the source is shown {green 
solid) with the la errorbars [green dashed). Also plotted are the X and Y residuals after 
subtracting off the best fit velocity (middle) and the best fit acceleration curve [bottom). 
The X (East-West) and Y (North-South) position plots (top) have a (y-axis) range of 0'.'16 
and residual plots [middle, bottom) have a (y-axis) range of ±8 mas. 
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Fig. 3. — Measured positions and residuals as a function of time for IRS 16NW, a source 
that has an acceleration consistent with zero, but significantly below the maximum possible 
acceleration. See the caption in Figure [2] for more information. 
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Fig. 4. — Histograms of a^jaap showing the significance of the acceleration measurements in 
both the radial (to'p) and tangential {bottom) directions. While, the distributions show an 
offset from zero indicating a possible bias due to systematic errors, such as residual distortion, 
that are not well characterized, it appears that any biases are limited to the ~1(7 level. The 
only star with significant negative radial acceleration (>4(7) is SO-15 and it is assumed to be 
a real acceleration due to the gravity of the supermassive black hole. 
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Fig. 5. — Positions and proper motion vectors of the young stars in our sample. Candidate 
disk members are shown in red and non-disk members are shown in blue over-plotted on an 
LGS AO image in grey-scale. The names of the stars in the primary sample are shown in 
the left panel and the complete extended sample is shown in a zoomed-out view in the right 
panel. The position of Sgr A* is marked with a black cross. 
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Fig. 6. — The significance of observed limits for the plane-of-the-sky acceleration. For a given 
projected radius, there is a maximum allowed acceleration (dashed line). If the measured 
accelerations from polynomial fitting, shown as So" upper limits on the y-axis, are less than 
the maximum allowed acceleration (below the dashed line), then significant constraints can 
be placed on the line-of-sight distance, z, and subsequently the orbital parameters of the 
star. SO-15 has a significant detection of non-zero acceleration and is plotted with its la 
errorbars. 
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Fig. 7. — The range of allowed orbital parameters for IRS 16SW as determined from the 
observed two-dimensional position in the plane of the sky, the three-dimensional velocity, 
and the acceleration. The probability distribution for each orbital parameter is determined 
by sampling from a gaussian distribution for each of the observed quantities and analytically 
converting to the standard orbital elements. High density (dark) regions represent the most 
probable values for each orbital parameter and the resulting la and 2a contours are shown 
as black lines. 
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Fig. 8. — The range of allowed eccentricities (e), inclinations (i), and angles to ascending 
nodes (Q) as determined by our orbit analysis for three example stars. The range of z values 
(horizontal axis) extends to all possible bound orbits for the star. The probability density 
function is shown in color with the la and 2cr contours drawn shown as black lines. SO-15 has 
a measured acceleration that is significantly different from zero. IRS 16C has an acceleration 
upper limit that is less than the maximum allowed acceleration, and thus a lower limit on 
the line-of-sight distance, \z\. IRS 16CC has no significant acceleration limit, but has a high 
velocity that is always larger than the circular velocity, thus prohibiting circular orbits. Also, 
by assuming the star is bound, the direction of the normal vector to IRS IGCC's orbital plane 
is restricted to a low inclination. 
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Fig. 9. — The resulting distribution of orbital phases for all stars when assuming either a 
uniform acceleration prior (black) or a uniform z prior (gray) and then imposing the mea- 
sured accelerations. The uniform z prior shows a strong bias towards an orbital phase of 0, 
which corresponds to periapse; while the uniform acceleration prior shows a more uniform 
distribution. 
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Fig. 10. — The orientation of three stars' orbital planes as described by the probability 
distribution of the planes' normal vector projected onto the sky as viewed from Sgr A*. 
Colors indicate the probability density for a star's normal vector to point at each pixel on 
the sky. The constraint on the stars' normal vectors are set by {top: SO-15) a measured 
acceleration; [middle: IRS 16C) a significant acceleration limit; [bottom: IRS 16CC) the 
star's high velocity and assuming the orbit is bound. 
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Fig. 11. — The la contours of all stars' probability distribution functions for the orientation 
of their orbital planes. This shows the distribution of stellar orbit orientations around the 
sky. The primary sample is plotted on the top and if there are degenerate solutions for a 
given star, then one solution is plotted with a solid line and the other with a dashed line. 
Additional sources found only in the secondary sample are plotted on the bottom and are 
plotted with dashed lines as there are no acceleration constraints and each star has a single 
solution with large uncertainties. We note that the orientation of the projection shown 
in this figure is rotated by 180° with respec t to that shown in earlier publications (e.g. 



Eisenhauer et al.ll2005l : iNayakshin et al.ll2006l ) in order to more easily see the region around 



the proposed disks. 
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Fig. 12. — The density of normal vectors to the orbital planes of the stars in our primary 
(top) and extended {bottom) samples. Densities are indicated in colors (stars deg~^) on a 
linear scale and the peak indicates an over-density of stars with similar orbital planes. Over- 
plot ted in black a r e the candidate orbital planes as proposed bvlLevin fc Beloborodovl (120031 ) 
and iGenzel et al.l (120031 ) with updated values from iPaumard et al.l (120061 ) for the candidate 
plane normal vector and uncertainties (solid black) and the disk thickness (dashed black) 
shown as solid angles of 0.05 sr and 0.09 sr for the clockwise and counter-clockwise disks 
respectively. 
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Fig. 13. — Positions for all candidate disk members in the disk plane from our primary 
(left) and extended (right) samples. The field of view for this study is projected onto the 
disk plane and shows the outer (dash) and inner (dot-dash) boundaries. For each disk 
candidate, a contour shows the star's position within the disk for all orbital solutions within 
10° of the disk plane. The color scale shows the probability density function for each star's 
position in the disk, normalized by the likelihood of disk membership. This normalization 
shows stars with a higher and lower likelihood of disk membership as darker red or lighter 
yellow, respectively. 
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Fig. 14. — The distribution of n-uncertainties as expressed by the area of the la region in 
which a star's normal vector can point. The uncertainties are shown both for the sample in 
this work [black) and for the stars w ith only three-diniensio nal position and two-dimensional 
velocity information extracted from iPaumard et al.l (120061 ) that are used in the search for a 
second disk (gray). 
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Fig. 15. — The radial distribution of stars within the disk plane for the extended sample. 
The best fit line is shown (dashed) and was constructed by excluding the first data point and 
the last three data points where field of view limitations may affect the distribution. 
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Fig. 16. — The distribution of eccentricity lower limits as determined from individual stellar 
orbits, excluding SO- 14. The top panel shows the 99.7% confidence lower limit from all pos- 
sible orbital solutions for candidate disk members {red circles) and non-disk members [blue 
squares). Stars from the primary sample {filled) and stars added in the extended sample 
(unfilled) are both shown. Sources in only the extended sample have less constrained eccen- 
tricities due to their larger velocity uncertainties. The bottom panel shows the candidate 
disk members 99.7% confidence lower limits after restricting the orbital solutions to those 
with normal vectors within 10° of the disk. By assuming disk membership, the range of 
eccentricities is more restricted for the candidate disk members. 
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Fig. 17. — Combined probability distribution for the candidate disk stars' eccentricity vec- 
tors. The eccentricity vectors for orbital solutions with normal vectors within 10° of the 
disk's normal vector are projected onto the disk plane. The la and 2a confidence-level 
contours are shown in black. 
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Fig. 18. — The improvement in positional accuracy at large radii as a result of correcting 
geometric distortion in speckle data sets. To characterize the systematic positional uncer- 
tainty, we take each star at each epoch and calculate the residual positional offset, which is 
defined as the difference between the measured position and the position as determined by 
the best fit velocity (x = + v * At) . Then the RMS of the residuals is calculated across 
all epochs for each star. All stars' resulting RMS values are sorted by the distance between 
the star and Sgr A* (which was at the center of the images) and then averaged over radius 
bins of 0'.'3. The radial trend is shown for data prior to the new distortion correction (top) 
and after the new distortion correction (bottom). 
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Fig. 19. — Map of the positional differences between stars observed near-simultaneously with 
NIRC and NIRC2. The maps are plotted in the original NIRC detector coordinates and show 
residuals before (left; a) and after (right; b) the NIRC-reimager distortion solution. 
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Fig. 20. — Distribution of the residuals before (gray) and after (black) correcting for the 
NIRC image converter distortion. Residuals are calculated by comparing a star's position 
in each NIRC image stack to the position in the LGS A0/NIRC2 image. These residuals 
are further reduced in the final image because the stacks are dithered small amounts on the 
detector and residual distortion can be averaged out if it is randomly oriented over the scale 
of the dither. 



